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The  new  matrix  will  represent  the  same  structure  as  the  original  matrix.  This  operation, 
which  amounts  to  the  partitioning  of  the  original  stiffness  matrix,  is  equivalent  to  the  elimina¬ 
tion  process  used  to  reduce  the  number  of  unknowns  in  a  set  of  simultaneous  equations  such  as 
Eq.  (2),  and  is  referred  to  as  "Elimination."  The  new  matrix  is  called  an  "Effective  Matrix." 

To  illustrate  the  elimination  process  consider  the  force-displacement  equation  for  unit  01 
in  Fig.  2; 


- 

- 

^0101 

^0101 

^0102 

^0102 

^0103 

^0103 

F 

^0104 

^0104 

^0105 

= 

K 

^0105 

^0106 

^0106 

^0304 

^0304 

F 

^030S 

^0305 

^0306 

^0306 

(4) 


where  F  and  U  represent  three  components  of  joint  force  and  joint  displacement,  respectively. 
Displacements  Uq^q^,  Uqjqz  ^0103  eliminated  from  the  equations  (these  are  at  the 

joints  not  common  to  other  units)  and  Eq.  (4)  is  reduced  to 


- 

• 

F’ 

^0104 

^0104 

F’ 

^  0105 

^0105 

0106 

K® 

^0106 

F‘ 

^0304 

i\ 

^0304 

^0305 

^0305 

^0306 

^0306 

The  original  matrix  (K)  of  size  27  x  27  is  thus  reduced  to  the  effective  matrix  (K®)  of  size  18  x  18. 
During  this  process,  the  force  vector  F  is  also  modified  to  F'.  The  details  of  this  modification 
are  described  in  Chapter  4. 

A  similar  effective  matrix  can  be  formed  for  unit  2  and  superimposed  on  that  for  unit  1.  The 
resulting  matrix  represents  the  structure  shown  in  Fig.  5.  At  this  point,  the  fact  that  displace¬ 
ments  and  must  equal  zero  is  taken  into  account  by  discarding  the  columns  and  rows 

in  the  matrix  corresponding  to  these  displacements.  The  matrix  for  the  two-unit  combination  is 
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then  25  x  25.  The  displacements  at  joints  0105  and  0104  and  the  remaining  displacement  at  0106 
are  now  eliminated,  leaving  an  effective  18  x  18  matrix. 

Finally,  the  effective  matrix  for  unit  3  is  superimposed.  This  does  not  increase  the  size  of 
the  matrix  since  no  new  joints  are  introduced.  Joints  0204,  0205  and  0206  are  now  eliminated 
leaving  a  9  x  9  matrix.  Displacements  and  yQ3Q5  are  all  zero  and  the  corresponding 

terms  in  the  matrix  may  be  discarded.  This  leaves  only  a  5  x  5  matrix,  which  is  finally  solved 
for  the  displacements  at  joint  0304  and  the  remaining  displacements  at  joint  0305. 

By  the  simple  process  of  substitution,  all  other  displacements  may  be  determined.  Having 
displacements,  all  bar  stresses  and  reactions  are  easily  obtained. 

In  comparison  to  the  method  described  above,  a  direct  solution  without  subdivision  of  the 
structure  would  have  required  the  solution  of  48  simultaneous  equations.  The  great  advantage 
of  the  method  is  now  obvious  when  applied  to  large  structures  having  many  joints. 


C.  INPUT  DATA 

It  is  desirable  to  prepare  the  input  data  on  format  input  data  sheets  before  punching  them 
onto  cards.  The  restrictions  of  the  STAIR  ANALYSIS  PROGRAM  and  the  items  of  data  required 
for  each  step  of  a  problem  are  listed  and  discussed  in  the  following  sections. 

Two  decks  of  cards  have  to  be  prepared  for  the  STAIR  System:  a  structural  data  deck  and 
a  load  data  deck.  These  decks  are  handled  differently  during  the  running  of  the  problem. 

The  structural  data  deck  is  input  directly  to  the  STAIR  System.  This  deck  begins  with  a 
block  of  cards  punched  from  an  "Initial  Data  Sheet."  Subsequent  blocks  of  cards  correspond  to 
consecutive  steps  in  the  ANALYSIS  PROGRAM;  the  first  card  of  each  block  calls  the  routine 
which  executes  the  step,  and  the  following  cards  provide  the  data  needed  for  the  execution. 

The  load  data  deck  is  input  to  the  STAIR  Load  Program  (SLOP).  This  program,  which  must 
be  run  before  the  ANALYSIS  PROGRAM,  receives  and/or  computes  concentrated  joint  loads  for 
various  types  of  loading  on  the  structure.  The  loads  are  assembled  into  groups  of  loads  applied 
to  individual  units  of  the  structure.  SLOP  is  stored  on  tape  for  use  by  the  STAIR  System. 

Each  item  of  data  entered  on  the  format  sheets  is  either  an  integer  number  (such  as  a  quan¬ 
tity  of  joints  or  bars)  or  a  fixed  number  (e.g.,  a  joint  coordinate  or  a  bar  area).  The  format 
sheets  prescribe  a  field  for  each  item  of  data.  A  field  is  the  maximum  number  of  characters 
an  item  may  contain. 

Integer  numbers  are  always  entered  at  the  right  of  their  field  and  are  never  written  with  a 
decimal  point.  Fixed  numbers  may  be  entered  in  one  of  two  ways:  with  or  without  a  decimal 
point.  If  the  decimal  point  is  specified  as  a  character  in  the  item  of  data,  the  number  may  be 
written  anywhere  within  its  field.  If  the  decimal  point  is  omitted,  the  number  must  be  positioned 
about  the  dashed  line  in  its  field,  with  the  fractional  part  to  the  right  of  the  dashed  line.  The  lo¬ 
cation  of  dashed  lines  is  written  into  the  machine  program  as  indicated  on  the  data  sheets  con¬ 
tained  herein.  If  there  is  no  dashed  line  In  the  field  of  a  fixed  number,  it  may  be  assumed  at 
the  right  of  the  field. 
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F.  INPUT  DATA  FOR  THE  STAIR  SYSTEM 


t 

c 


( 


( 


c 


The  input  parameters  required  by  each  step  of  the  STAIR  operations  are  described  below. 
The  restrictions  on  the  magnitude  of  these  parameters  and  operations  are  given  where  appli¬ 
cable.  The  formats  in  which  these  data  parameters  are  prepared  will  be  illustrated  in  the  fol¬ 
lowing  data  sheets.  Data  sheets  for  each  step  are  arranged  in  the  same  order  as  the  example 
for  the  ANALYSIS  PROGRAM.  The  same  example  structure  of  Fig.  2  is  used. 

Initial  Input  (refer  to  data  sheet,  p.  44) 

(a)  Problem  identification  number  (IDP). 

(b)  Number  of  units  into  which  the  given  structure  is  divided  (NUT) . 

(c)  Number  of  different  loading  conditions  to  which  this  given  structure 
is  subjected  (NK). 

(d)  Number  of  reaction  joints  or  supports  in  structure. 

(e)  Elastic  modulus  of  bar  material  E  kips  per  sq.  in. 

(f )  List  of  reaction  joint  numbers  and  the  coordinate  directions  in 
which  displacement  is  restricted  (1 -immovable;  0-free). 

Size  Restrictions 

(a)  Number  of  units  <  99. 

(b)  Number  of  reaction  joints  <  125. 

(c)  Loading  conditions  <  28. 

STRIX  (refer  to  data  sheet,  p.  45) 

(a)  Step  number  in  ANALYSIS  PROGRAM.  It  precedes  every  operation 
code,  as  the  001  on  p.  45  precedes  01  STRIX,  and  the  005  precedes 
04  MADD  on  p.  49.  It  is  the  first  card  of  a  card  block  calling  the 
routine  which  executes  the  step  as  mentioned  in  the  section  on  input 
data.  (Step  numbers  must  start  at  01  and  proceed  sequentially  to 
the  end.) 

(b)  Identification  number  of  the  unit  (NUU)  to  which  this  STRIX  operation 
applies.  On  p.  45  NUU  is  01;  later,  units  02  and  03  will  be  formed 
by  MATRO.  In  case  some  of  the  units  are  not  identical,  NUU  is  used 
to  identify  the  different  STRIX  operations  for  different  units  when 
STRIX  is  called  more  than  once. 

(c)  Number  of  joints  in  the  unit  (NJTU) . 

(d)  Number  of  bars  in  the  unit  (NMU). 

(e)  Joint  sequence.  This  serves  as  a  reminder  only,  and  is  not  a  part 
of  the  input  data. 

(f)  Joint  cards  giving  the  joint  numbers  [see  Operation  Restriction  (b)] 
contained  in  the  unit  in  sequence*  and  the  coordinates  of  the  joints 
in  feet  (reproduced  from  SLOP). 

(g)  Bar  cards  are  reproduced  from  SLOP.  Bars  which  are  common  to 
two  units  are  listed  separately  for  this  input  (see  p.  46). 


♦The  sequence  can  be  of  any  convenient  form,  but  the  joints  that  are  to  beveliminated  by  EFFRIX 
must  be  placed  first  in  the  sequence. 
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size  Restrictions 


If  N  =  3x  (number  of  joints  in  the  unit),  NK  =  number  of  loading  conditions, 

(a)  t  1)  +  N  X  NK  <  16000  (see  table). 


NJTU 

NK 

59 

1 

58 

4 

57 

7 

56 

10 

55 

13 

54 

17 

53 

20 

52 

24 

51 

27 

(b)  Number  of  bars  <  600. 

(c)  Unit  number  <  total  number  of  units  in  the  structure. 

Operation  Restrictions 

(a)  In  no  case  may  all  the  bars  meeting  at  a  joint  lie  in  one  plane  unless 
the  joint  is  a  reaction  joint.  If  a  planar  joint  exists,  a  new  bar  out 
of  the  plane  (an  imaginary  member)  may  be  inserted  which  will  carry 
no  stress  but  will  allow  the  analysis  to  proceed.*  The  new  bar  may 
connect  the  planar  joint  with  any  other  joint  not  lying  in  the  plane. 

(b)  If  an  EFFRIX  operation  immediately  follows  STRIX  in  the  sequence 
of  the  ANALYSIS  PROGRAM,  those  joints  which  are  going  to  be  elim¬ 
inated  have  to  be  placed  first  in  the  row  of  joint  cards.  In  the  ex- 
ample,  joints  0101,  0102  and  0103  are  eliminated  by  the  following 
EFFRIX. 


EFFRIX  (refer  to  p.  47) 

(a)  Step  number  in  ANALYSIS  PROGRAM. 

(b)  Number  of  joints  (NE)  to  be  eliminated  from  the  row  sequence 
established  by  the  previous  operation  (STRIX  or  MADD). 

Size  Restrictions 


NE  =  number  of  joints  to  be  eliminated. 

NJ  =  number  of  joints  in  original  row  sequence  (number 
of  joints  in  the  unit) . 

NR  =  number  of  reaction  joints  among  joints  being  eliminated. 

NC  =  number  of  reaction  components  acting  at  the  eliminated 
reaction  joints. 


*  See  p.  41  for  a  more  detailed  explanation  of  a  planar  joint. 
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G.  ERROR  HALTS  IN  THE  BSS  LOADER 


The  STAIR  system  uses  the  old  BSS  loader  (commonly  called  the  10-card  loader)  and  not  the 
newer  BSS2  loader  because  of  incompatible  timing  and  binary  patches  that  are  made  to  the  loader 
by  the  MISLAM  routine  at  execution  time.  The  BSS  loader  has  known  errors  in  it;  however,  these 
do  not  affect  the  STAIR  system.  Because  this  loader  is  essentially  obsolete,  information  on  its 
operation  and  error  stops  is  difficult  to  obtain.  For  this  reason,  the  following  table  (as  contained 
in  an  early  IBM  709  FORTRAN  Operators  Manual  and  Lincoln  Memorandum  No.  2C-0318)  is  re¬ 
produced  here  for  the  users  convenience. 


Halt  Location 

(Octal)  Reason  for  Halt 


77525  Check  sum  error  on  cards. 

77755  Instructions  overlap  the 

Symbol  Table  of  the 
Loader. 

77400  Instructions  and  Data 

overlap. 

77731  More  than  20  subroutines 

are  missing. 

77437  EOF  at  the  card  reader. 

This  condition  does  not 
produce  a  Halt.  The 
Program  Light  is  not 
turned  on  but  the  Read 
Select  Light  is  turned  on. 
The  location  counter  will 
contain  77437. 

77747  Missing  subroutines. 


Procedure 

Press  START  to  accept  information. 

Terminate  loading.  Combination  of 
program  and  transfer  vectors  too  long. 
Rewrite  program. 

Terminate  loading.  Combination  of  in¬ 
structions  and  data  too  long.  Rewrite 
program. 

If  missing  subroutines  are  immediately 
available,  ready  in  the  card  reader, 
press  START  until  a  stop  at  77747  occurs. 

Press  START  to  read  cards. 


This  stop  indicates  the  TRANSFER  CARD 
has  been  reached.  It  is  caused  by  one  of 
the  two  occurrences  listed  below: 

(1)  Loading  has  been  completed,  but  at 
least  one  of  the  subroutines  called  for  is 
missing.  Location  77400  contains  the 
BCD  name  of  the  first  missing  subroutine, 
location  77401  the  second,  etc.  If  the 
missing  subroutine(s)  is  immediately 
available,  it  may  be  loaded  without 
starting  the  entire  loading  process  over 
again.  Place  another  TRANSFER  CARD 
(9  punch  in  column  1)  at  the  end  of  the 
routine(s),  ready  in  the  card  reader, 

and  press  START. 

(2)  The  TRANSFER  CARD  encountered 
is  really  a  premature  one  that  simply 
has  not  been  withdrawn.  Be  certain  that 
a  TRANSFER  CARD  is  the  last  card  at 
the  end  of  the  deck  and  press  START. 
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D.  DEFINITION  OF  COMMON  VARIABLES 


Common  storage  contains  the  following  list  of  variables  during  execution  of  aU  the  operating 
routines; 


IDP 

NSTEP 

E 

NREF 

NJREAC 

NUT 

NK 

NU 

NJT 

NJSEC 

NRJSEC 

NRCSEC 

JREAC 


NT  A  PE 

NROW 

STRDC 


problem  identification  number 

consecutive  number  of  steps  in  ANALYSIS  PROGRAM  being 
executed 

elastic  modulus  of  bar  material 

number  of  residual  and  reaction  matrices  written  on  tape  A6 
number  of  reaction  joints  in  structure 
number  of  units  in  structure 
number  of  loading  conditions 

number  of  unit  whose  matrix  is  in  primary  storage  area 
number  of  joints  in  primary  row  sequence 

number  of  joints  being  eliminated  from  primary  row  sequence, 
or  number  of  joints  in  secondary  row  sequence 

number  of  reaction  joints  among  joints  being  eliminated 
number  of  reaction  components  among  joints  being  eliminated 

list  of  all  reactions  acting  on  structure  (four  words  per  reaction 
joint  giving  joint  number  and  code  indicating  whether  or  not  com¬ 
ponents  of  reaction  act  in  x-,  y-  or  z-directions);  (0)  means  no 
reaction,  (1)  means  a  reaction 

tape  reading  error  indicator  (set  by  FIXTPE  and  sensed  by 
routine  calling  FIXTPE) 

list  of  joint  numbers  associated  with  matrix  in  primary  storage 
area 

primary  matrix  storage  area 


The  following  arrays  are  included  in  common  storage  by  all  routines  except  STRDC  (to  be  dis¬ 
tinguished  from  array  STRDC),  BARSOL  and  CHECK, 


NROSEC  secondary  row  sequence,  associated  with  matrix  in  secondary 
storage  area 

LRSEC  list  of  reactions  among  joints  being  eliminated 
STSEC  secondary  matrix  storage  area 

Additional  parameters  and  arrays  may  be  included  in  common  storage  by  the  use  of  certain  op¬ 
erating  routines,  but  these  parameters  do  not  affect  the  continuity  of  operation  among  routines. 


E.  TAPE  STORAGE  FORMATS 

Storage  formats  for  the  input  and  buffer  tapes  attached  to  DSC  A  are  given  below.  The 
format  of  the  program  tape  B4  will  be  described  with  MISLAM. 

Tape  A6  —  Residuals  and  reaction  matrices  in  two  binary-record  groups. 

Record  1  -  IDOUT,  NJT,  NJSEC,  NRJSEC,  NRCSEC 
Record  2  -  NROW,  LRSEC,  STSEC 

Tape  A7  -  Matrix  equation  storage  in  three  binary-record  groups. 

Record  1  -  NHED,  NJT,  NJSEC,  NRJSEC,  NRCSEC,  NU 

Record  2  -  NROW.  STRDC 

Record  3  -  NROSEC,  LRSEC,  STSEC 


Tape  A8  —  Load  data  with  consecutive  binary  records  containing  the  load 

matrix  for  consecutive  units  of  the  structure.  The  first  word  of 
each  record  is  a  heading  NHED. 

Tape  A9  —  Structural  data  in  groups  of  BCD  records. 

Following  an  initial  group  of  two  data  records,  the  records  for 
the  step  number  and  the  operation  number  follow,  each  of  which 
is  followed  by  the  data  records  for  the  step. 

Records  1  and  2  —  General  structure  data 
Record  3  —  Step  number,  operation  number 
Record  4. ...  n  -  Data  for  this  step 
Record  n+1  —  Step  number,  operation  number 

NHED  identifies  a  record  group  u  as  32767-u. 

IDOUT  identifies  a  two-record  group  as  a  residual  (0)  or  reaction  matrix  (1). 

F.  STAIR  OPERATING  ROUTINES 

The  nine  operating  routines  are  described  on  the  following  pages.  Each  description  has  the 
following  format: 

(a)  Operation  —  definition  of  the  matrix  operation  performed. 

(b)  Mathematical  method  —  brief  description  of  technique  for  performing 
the  operation. 

(c)  Initial  condition  of  core  —  common  variables  which  must  be  present  to 
maintain  continuity. 

(d)  Error  stops. 

(e)  Output  -  to  tape  only;  on-line  output  consists  of  error  diagnostics  and 
tracing  information. 

(f )  Final  condition  of  core  —  common  variables  which  are  changed,  affecting 
continuity. 

(g)  Addressing  diagrams. 

(h)  Block  diagrams  —  flow  charts. 


1.  STRDC 

Operation:  Forms  a  stiffness  matrix  equation  in  core  from  input  defining  the  geometry  of 
a  structure  and  its  load  matrix. 

Mathematical  Method:  Terms  of  the  matrices  for  single  bars  are  computed  and  superim¬ 
posed  on  the  array  of  the  final  matrix.  The  location  of  the  final  matrix  array  of  the  bar  con¬ 
necting  joints  a  and  b  is  determined  by  the  positions  of  a  and  b  in  the  list  of  joint  numbers 
NROW.  The  load  matrix  for  the  structure  is  read  directly  from  tape  A8,  and  is  stored  imme¬ 
diately  after  the  stiffness  matrix.  A  vector  check indicates  if  all  the  bars  meeting  at  any  joint 
lie  in  a  plane;  such  a  geometry  would  result  in  a  singular  stiffness  matrix. 

Initial  Condition  of  Core:  Common  variables  E,  NUT  and  NK  must  be  predetermined. 


1‘See  footnote  on  p.41  for  an  explanation  of  this  check. 
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EFFRIX  ADDRESSING  DIAGRAMS 
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<,0 


NRA  =  3  NJSEC  -  NRCSEC 
NRB  =  3  (NJT-  NJSEC) 

SIZE  RESTRICTIONS: 

NRA  (NRA  +  l)/2  +  NRA  (NRB  +  NK) «  9500 

NRA  (NRA  +  1)<  3  NJSEC  (3  NJSEC  +  lyz  +  9  NJSEC  (NJT  -  NJSEC) 
(3  NJT  +  NRB  +  I )  3  NJSEC/2  >  (NRB  +  1 )  NRB/2 
NRA>0 
NRJSEC4  10 


MATRIX  RESIDUAL  STORAGE  BY  EFFRIX 
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Experience  with  the  STAIR  program  over  the  past  few  months  has  disclosed  a  number  of  errors 
in  the  printed  manual.  Rather  than  list  the  errors  on  a  sheet  and  rely  on  each  holder  of  the 
manual  to  make  the  appropriate  corrections,  corrected  sheets  for  insertion  into  the  manual  as 
follows  are  herein  provided: 


Page  15; 

Page  16: 

Page  17; 

Page  21 : 

Page  25; 
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Page  80: 

Page  84: 
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Page  123; 


Problem  200,  SLOP  Input  Data  Sheet  'B,'  Structural  Data  I:  coordinate 
listing  corrected.  Remove  original  page  15. 

Problem  200,  SLOP  Input  Data  Sheet  'B,'  Structural  Data  II:  area 
listing  corrected.  Remove  original  page  16. 

Problem  200,  SLOP  Input  Data  Sheet  'B,'  Structural  Data  II:  area 
listing  corrected.  Line  out  page  17;  save  page  18. 

Problem  200,  SLOP  Input  Data  Sheet  'C,'  Loading  Data,  Case  #3: 
loading  corrected.  Line  out  page  21;  save  page  22. 

Printout:  print  error  corrected.  Remove  original  page  25. 

Printout;  print  error  corrected.  Remove  original  page  26. 

Program  description:  equations  corrected.  Line  out  page  37; 
save  page  38. 

Program  description;  equations  corrected.  Line  out  page  39;  save 
page  40. 

Program  description:  equations  corrected.  Line  out  page  80;  save 
page  79. 

Program  block  diagram;  equations  corrected.  Line  out  page  84; 
save  page  83. 

DEFSOL  Matrix  Addressing  Diagram:  equations  corrected.  Line 
out  page  118;  save  page  117. 

DEFSOL  Matrix  Storage  Diagram,  Matrix  Storage  Area;  equations 
corrected.  Line  out  page  119;  save  page  120. 

Program  block  diagram:  equations  corrected.  Line  out  page  123; 
save  page  124. 
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20UNES 
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30 

0201  0204- 

30 

0202  0203 

30 

0202  0204- 

30 

0207  0205 

30 

0202  020£> 

30 

0203  0204. 

30 

0203  0205 

30 

0204  0205 

30 

0204  020& 

30 

0204  0301 

30 

0204  0302 

30 

0204  0303 

30 

0205  0206 

30 
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NK  =  number  of  different  loading  conditions. 

A  =  3NE  -  NC. 

B  =  3(NJ-NE). 

C  =  3NE. 

(a)  NR  <  10. 

(b)  +  a(B  +  NK)  <  9500. 

(c)  A{A  +  1)  <  (C  X  5_+J.  +  c  X  B], 

(d)  <  [C  X  +  C  X  B). 

(e)  A  >  0 

MATOUT  (refer  to  p.  47) 

(a)  Step  number. 

(b)  Tape  location  (Loc.)  In  which  a  matrix  is  to  be  stored.  Matrices 
are  stored  on  tape  in  bloc)cs.  These  are  numbered  consecutively, 
the  tape  location  being  the  number  given  to  the  block. 

Operation  Restriction 

Only  one  matrix  may  be  stored  in  a  tape  location. 

MATRO  (refer  to  p.  48) 


(a)  Step  number. 

(b)  Identification  niunber  of  new  unit  (NUN). 

(c)  Angle  through  which  old  unit*  must  be  rotated  in  the  x-y  plane  to 
arrive  at  the  position  of  the  new  unit.  Angle  (EPS)  measured  in 
degrees  and  positive  when  clockwise,  as  viewed  by  an  observer 
on  the  positive  z-axis.  In  the  STAIR  system,  the  rotation  must 
be  about  the  z-axis. 

(d)  Row  sequence  of  the  resulting  matrix  after  rotation  [NROW  (1) 
to  NROW  (NJT)J. 

(e)  If  the  matrix  being  rotated  is  a  reduced  matrix,  i.e.,  some  joints 
have  been  eliminated  from  it  by  EFFRIX  in  a  previous  step,  the 
inverse  matrix  corresponding  to  these  eliminated  joints  is  stored 
in  the  secondary  matrix  storage  —  then,  not  otherwise.  The  new 
row  sequence  of  the  eliminated  joints,  i.e.,  the  new  row  sequence 
for  the  new  inverse  matrix,  is  NROSEC  (1)  to  NROSEC  (NE). 

Size  Restriction 

(a)  New  unit  number  <  number  of  units  in  the  structure. 

Operation  Restriction 

(a)  Matrices  of  identical  units  from  which  joints  )iave  been  eliminated 
may  be  formed  by  MATRO  except  when  one  of  the  eliminated  joints 
of  the  old  unit  has  a  corresponding  joint  in  the  new  unit  which  is  a 
reaction  joint,  or  vice-versa. 

(b)  MATRO  can  not  follow  a  double  EFFRIX  sequence. 


{ 


♦Matrix  of  unit  in  core. 
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DEFSOL  (refer  to  p.  “'•si) 


(a)  Step  number. 

(b)  List  of  the  number  of  joints  which  are  prefixed  by  each  unit  number 
(lUAOR).  For  example,  there  are  9  joints  in  unit  01,  but  only  6 
joints  are  prefixed  by  unit  number  01.  (Notes  regarding  thll  In 
SLOP  apply  here  also.) 


Size  Restrictions 

NJ  =  number  of  joints  in  final  row  sequence. 

NR  =  number  of  reaction  joints  among  final  joints. 

NC  =  number  of  reaction  components  among  final  joints. 
NK  =  number  of  loading  conditions. 

NJT  =  total  number  of  joints  in  structure. 

NRT  =  total  number  of  reaction  joints  in  structure. 

A  =  3NJ  -  NC. 

(a)  NR  <  10. 

(b)  A  X  (A  +  l)/2  +  A  X  NK  <  9500. 

(c)  A(A  +  1)  <  16,000. 

(d)  4  X  NK(NJT  +  NRT)  <  16,000. 


Operation  Restriction 


DEFSOL  may  be  called  only  when  the  final  matrix  of  the  structure  is  stored  in  the  computer 
core. 


Output 


(a)  Joint  displacement  (in  Inches). 

(b)  Reactions  (in  kips). 


BARSOL  (refer  to  p.  54) 

(a)  Step  number. 

(b)  Number  of  bars  whose  stresses  are  to  be  computed  (NM). 

(c)  Number  of  joints  to  which  the  bars  listed  in  (b)  are  connected  (NJ). 

(d)  List  of  bar  numbers  of  the  bars  whose  stresses  are  to  be  computed 
(reproduced  from  SLOP). 

Note:  STAIR  will  compute  any  desired  number  of  bar  stresses  by 
repeated  use  of  BARSOL. 

(e)  Joint  cards,  giving  the  joint  numbers  and  coordinates.  Only  those 
joints  which  appear  in  the  list  of  bar  numbers  (p.  54)  need  be  in¬ 
cluded  (reproduced  from  SLOP). 

Size  Restrictions 

NK  =  number  of  loading  conditions. 

(a)  NJ<  800. 

(b)  NM  <  1000. 

(c)  NK  X  NM  <  4300. 
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SIZE  RESTRICTIONS; 
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OEFSOL  MATRIX  ADDRESSING  DIAGRAM 


REAaiON  EQUATIONS  (STSEC) 


SIZE  LIMITATIONS; 
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DEFSOL  MATRIX  STORAGE  DIAGRAM 


MATRIX  STORAGE  AREA 
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CHAPTER  1 
INTRODUCTION 


An  analysis  of  a  structure  is  a  determination  of  its  deflections,  reactions  and  internal 
stresses  under  load.  For  purposes  of  STAIR,  a  structure  is  defined  as  a  three-dimensional  sys¬ 
tem  of  prismatic  bars,  connected  at  frictionless  pin  joints  and  supported  in  a  stable  position  by 
reactions.  The  only  external  loading  considered  is  a  system  of  concentrated  forces  applied  to 
the  joints.  When  the  number  of  joints  making  up  such  a  framework  becomes  large,  the  conven¬ 
tional  methods  of  analysis  become  too  cumbersome  for  slide  rule  or  desk  calculator  solutions, 
and  a  high-speed  machine  analysis  is  necessary.  If  it  is  necessary  to  use  electronic  digital  com¬ 
puters,  it  is  desirable  to  consider  less  common  methods  of  analysis,  which  may  be  more  suitable 
than  conventional  procedures  for  use  on  a  machine. 

The  stiffness  matrix  method  of  analysis,  which  is  the  basis  of  STAIR,  is  such  a  method.  The 
method  has  two  essential  steps:  (1)  the  forming  of  a  stiffness  matrix  for  the  structure  from  the 
geometric  properties  of  the  structure,  and  (2)  the  solving  of  a  matrix  equation,  which  is  a  set  of 
simultaneous  equations  formed  from  selected  coefficients  of  the  stiffness  matrix  and  the  applied 
loads.  The  set  of  equations  describes  the  equilibrium  of  forces  at  the  movable  joints  in  terms 
of  their  unknown  displacements.  The  solution  of  the  equations  is  a  set  of  displacements,  which 
can,  in  turn,  be  used  to  solve  for  the  bar  stresses  and  reactions  of  the  structure. 

The  characteristic  of  this  method  that  makes  it  amenable  to  machine  computation  is  its  re¬ 
petitive  nature.  The  form  of  the  stiffness  matrix  is  the  same  for  each  member  and  a  stiffness 
matrix  for  a  structure  is  merely  the  sum  of  the  matrices  for  the  individual  members.  Solving 
a  set  of  simultaneous  equations  is  usually  accomplished  by  means  of  some  sort  of  algorithm, 
which  is  also  a  repetitive  procedure.  A  second  advantage  of  this  method  is  its  independence  of 
the  degree  of  redundancy.  Additional  redundancies  do  not  require  more  machine  time.  Finally, 
the  procedure  yields  all  the  joint  displacements  as  the  first  available  result,  and,  thus,  is  es¬ 
pecially  suited  to  displacement  problems. 

The  one  characteristic  of  the  stiffness  matrix  method  of  analysis  that  is  somewhat  incon¬ 
venient  and  which  makes  it  unsatisfactory  for  use  in  other  than  machine  computation  is  the  large 
number  of  simultaneous  equations  that  must  be  handled.  In  general,  there  are  three  unknown 
displacements  at  each  joint  of  a  three-dimensional  structure;  therefore,  three  equations  must 
be  written  at  each  joint.  This  becomes  an  inconvenience,  even  in  machine  computation,  when 
larger  structures  with  more  joints  are  considered.  The  STAIR  procedure  has  been  developed  to 
break  down  large  structures  into  small  parts  because  the  simultaneous  equations  required  for 
large  structures  do  not  fit  into  the  available  memory  space  in  many  computers  and,  more  impor¬ 
tantly,  require  an  excessive  amount  of  computer  time  if  solved  by  direct  methods.  This  enables 
the  programmer  to  remain  within  the  size-time  restrictions  of  the  modern  computer. 

The  analysis  of  any  three-dimensional  framework  by  this  method  depends  on  the  proper  com¬ 
bination  of  nine  "Basic  Operations"  which  have  been  prepared  as  self-contained  computer  "sub¬ 
routines."  The  internal  functioning  or  derivation  of  these  operations  need  not  concern  the  struc¬ 
tural  analyst,  but  he  must  have  clearly  in  mind  their  functions  and  limitations. 

Among  the  requirements  of  each  subroutine  are  certain  numbers  which  describe  the  geometry 
of  the  structure  and  parameters  which  are  present  for  a  given  structure  to  control  the  manipula¬ 
tion  of  geometrical  data  by  the  subroutines. 
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The  ANALYSIS  PROGRAM  described  in  Chapter  2  is  the  device  for  effecting  an  analysis  by 
this  procedure.  It  is  essentially  a  list  of  the  names  of  the  necessary  subroutines  arranged  in 
proper  sequence  for  a  particular  structure.  Associated  with  each  subroutine  is  its  data  and 
control  parameters.  Rigid  rules  are  presented  for  each  subroutine  and  its  data.  If  these  are 
followed  exactly,  it  is  easy  to  analyze  almost  any  three-dimensional  pin-jointed  framework. 

In  the  following  chapters,  the  stiffness  matrix  method  of  analysis  is  explained  in  some  de¬ 
tail,  the  ANALYSIS  PROGRAM  is  described,  the  rules  for  use  of  the  subroutines  are  listed,  the 
possible  errors  and  pitfalls  are  listed  for  each  operation  and  a  sample  solution  is  presented. 
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CHAPTER  2 
USER'S  MANUAL 

A.  INTRODUCTION 

Various  techniques  exist  for  analyzingthe  structural  behavior  of  a  stable,  three-dimensional, 
pin-jointed  framework.  Generally,  statically  determinate  structures  are  analyzed  by  solving 
equations  of  equilibrium  at  joints  or  from  cut  sections,  whereas  indeterminate  structures  require 
the  solution  of  additional  simultaneous  equations  that  are  related  to  the  redundancies  of  such 
structures. 

These  methods  were  developed  for  hand  and  slide  rule  computation;  in  many  cases  the  re¬ 
sults  are  at  an  acceptable  level  of  precision.  However,  the  introduction  of  high-speed,  exter¬ 
nally  programmed,  electronic  computing  machines  has  opened  new  channels  for  analysis.  In 
particular,  it  has  permitted  the  engineer  to  use  rigorous  mathematical  models  of  structures, 
which  formerly  could  not  be  used  on  large  frameworks  due  to  the  sheer  volume  of  computations 
required. 

The  STAIR  system  has  been  developed  to  use  the  stiffness  matrix  technique,  which  handles 
a  structure  in  terms  of  the  stiffness  of  its  members  and  the  deflections  of  its  joints.  This  sys¬ 
tem  is  outlined  in  Sec.  B. 

The  geometry  of  a  structure  is  defined  as  the  position  and  orientation  of  its  joints  and  mem¬ 
bers  relative  to  some  fixed  coordinate  system,  the  cross-sectional  area  and  modulus  of  elasticity 
of  its  members  and  the  location  and  direction  of  its  reactions.  This  is  a  necessary  ingredient 
in  any  structural  analysis  technique. 

The  size  of  the  core  memory  of  the  newest  machines,  while  large,  is  still  the  limiting  factor 
in  the  size  of  structure  which  can  be  analyzed.  Since,  in  many  cases,  the  structure  under  con¬ 
sideration  may  be  several  times  larger  than  the  limits  of  the  memory,  STAIR  is  arranged  to 
permit  the  division  of  the  structure  into  small,  internally  stable  units,  which  can  be  partially 
analyzed  and  then  combined  for  a  final  solution.  This  ability  imposes  a  definite  nomenclature 
on  the  geometrical  arrangement  of  the  structure. 

The  nomenclature  used  in  STAIR  is  shown  in  Figs.  1  and  2.  If  the  structure  is  small  enougli, 
it  may  be  considered  as  a  single  unit  (Fig.  1)  and  analyzed  directly.  However,  as  the  number  of 
joints  (and  therefore  equations)  increases,  it  is  necessary  to  subdivide  the  structure,  as  shown 
in  Fig.  2.  It  should  be  noted  that  STAIR  is  particularly  useful  with  structures  showing  a  rotational 
symmetry  around  the  z-axis.  Units  are  numbered  consecutively  starting  with  01.  Joints  are 
given  a  four-digit  number,  the  first  two  being  the  unit  number  with  which  it  is  identified.  Thus, 
joint  0112  in  Fig.  1  is  the  12th  joint  in  unit  01,  and  joint  0206  in  Fig.  2  is  the  6th  joint  in  unit  02. 

To  avoid  duplication,  a  joint  is  given  only  one  number  even  though  it  may  appear  in  two  units. 
Thus,  in  Fig.  2,  joint  0104  is  shown  both  in  unit  01  and  unit  02.  Bars  are  designated  by  eight¬ 
digit  numbers,  which  consist  of  the  numbers  of  the  joints  which  the  bar  connects.  The  joint  with 
the  lower  number  is  listed  first  by  convention.  Thus,  the  bar  connecting  joints  0104  and  0105 
is  denoted  as  0104-0105.  In  this  case,  where  the  bar  is  common  to  two  units,  the  programming 
is  handled  by  using  one-half  of  the  cross-sectional  area  in  the  input  of  each  unit. 

The  procedure  for  solving  a  structure  with  STAIR  is  as  follows; 

1.  Load  data  are  prepared  for  as  many  cases  as  desired.  The  data,  including  dead  load 
information,  the  geometry  of  the  structure  and  the  various  live  load  conditions,  are  introduced 
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Into  the  STAIR  Load  Program  (SLOP);  the  output  of  this  program  la  a  data  tape  which  serves  as 
Input  to  the  STAIR  ANALYSIS  PROGRAM, 

2,  The  ANALYSIS  PROGRAM  Is  prepared,  listing  the  various  subroutines  required  for 
solution.  Geometrical  data  for  the  structure  are  again  Introduced,  Including  the  division  of  the 
structure  Into  smaller  units.  When  possible,  units  should  have  Identical  bar  arrangements  with 
rotational  geometric  symmetry,  since  this  will  reduce  Input  data  requirements  and  computational 
time.  Loads  and  reaction  joints  for  these  units  need  not  be  symmetric. 

3.  THE  ANALYSIS  PROGRAM,  the  geometrical  data  and  the  SLOP  output  tape  are  Intro¬ 
duced  to  the  computer,  which  will  form,  modify  and  solve  the  stiffness  matrices,  and  then  back- 
substitute  to  get  the  final  structural  solution.  The  results  of  a  complete  run  consist  of  the  joint 
displacements,  bar  stresses,  reactions  and  summation  of  internal  bar  forces  at  each  joint  for 
each  loading  condition. 

In  the  following  sections,  a  complete  program  for  the  example  structure  of  Fig.  2  is  set  up 
and  solved  on  the  IBM-7090.  Descriptions  of  the  input  data  and  program  arrangement  and  dis¬ 
cussion  of  the  size  restrictions  are  presented  at  each  step.  Although  the  structure  used  is  small 
enough  to  permit  solving  as  a  single  unit,  it  has  been  divided  into  rotationally  symmetric  units 
to  demonstrate  the  use  of  the  various  subroutines. 


B.  THE  STIFFNESS  MATRIX  METHOD 


The  stiffness  matrix  method  used  in  structural  analysis  is  the  solution  of  the  force  equilib¬ 
rium  equations  in  the  form  of  a  matrix  equation.  The  STAIR  System  is  mainly  a  matrix  opera¬ 
tion,  such  as  matrix  formation,  addition,  multiplication,  inversion,  elimination,  etc.  When  a 
final  stiffness  matrix  representing  the  structure  is  obtained,  it  is  solved  for  the  joint  displace¬ 
ments  (U)  corresponding  to  a  particular  loading  condition  (F). 

A  stiffness  matrix  (K)  is  a  numerical  representation  of  a  structure.  It  relates  the  joint  load 
(F)  and  joint  displacement  (U)  in  the  form  (F)  =  (K)(U).  Consider  the  bar,  0102-0104,  of  the 
example  structure.  Figure  3  shows  its  relation  to  the  xyz  coordinate  system  and  the  uvw  dis¬ 
placement  system.  Figure  4  shows  the  bar  after  it  has  undergone  positive  u,  v  and  w  displace¬ 
ments.  Note  that  the  subscript  1  refers  to  joint  0102,  and  subscript  2  refers  to  joint  0104. 

From  Fig.  4, 


f^l  =  =  F^a 

fyj  =  Fj  cos  <P2  = 
f^l  =  Fj  cos  <pj  =  F^y 


(1) 


where  ^xl'  ^1  and  fzl  are  the  x,  y  and  z  forces  at  joint  1  (0102),  f^^.  f^^  Y 

and  z  forces  at  joint  2  (0104), 


y  =  Tr 


are  the  direction  cosines  and  L  is  the  length  of  the  bar.  F^  can  now  be  related  to  the  uvw  dis¬ 
placements  in  the  following  manner; 

AE 

Fj  =  “17  [(«!  -  “2)  cos  <if>  j  +  (Vj  -  v^)  cos  If  2  +  -  w^)  cos  v>  j) 
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Fj  =  —  ((Uj  -  U2)o(  +  (v^  -  +  (Wj  -  w^jy)  , 

where  A  la  the  croaa-aectional  area  of  the  bar  and  E  is  its  modulus  of  elasticity. 

Substitution  into  Eq.  (1)  gives 

f^j  =  AE/Lfo^Uj  +  +  ayw^  -  a^U2  -  1 

fy^  =  AE/L(a^u^  +  +  ^yw^  -  a^U2  fiy'Nj] 

^zl  ""  AE/L(ayu^  +  ^yv^  +  y^w^  -  aya.^  -  -  y^W2] 

z  z 

fj^2  =  AE/L(-a  Uj  -  -  oyw^  +  a  U2  +  oi^w.^  +  ayw2] 

fy2  =  AE/L(-a^Uj  -  -  /3yw^  +  allu^  +  +  Pyvi^] 

fjj2  =  AE/L[-fl(yu^  -  pyv^  -  y^w^  +  ayU2  +  Pyy^  +  y^''>z^  •  (2) 

A  convenient  way  of  writing  Eq.  (2)  is  to  employ  matrix  notation  as  follows: 


or  (f)  =  (k)  (u). 

The  complete  representation  of  the  whole  structure  or  unit  can  be  obtained  by  superimposing 
(f)  =  (k)  •  (u)  for  all  bars,  giving  (F)  =  (K)  •  (U)  for  the  whole  structure  or  unit.  Then,  stiffness 
matrices  representing  adjacent  units  of  a  structure  may  be  superimposed  to  give  a  stiffness  ma¬ 
trix  representing  a  combined  structure. 

The  formation  of  a  stiffness  matrix  for  a  new  unit,  which  is  rotationally  symmetric  to  an¬ 
other  unit,  requires  only  a  transformation  of  the  stiffness  matrix  of  the  latter  unit.  The  trans¬ 
formation  in  the  present  STAIR  System  is  limited  to  rotation  of  the  plane  containing  two  coordi¬ 
nate  axes  (the  x  and  y)  about  the  third  (the  z  -axis) . 

The  size  of  the  stiffness  matrix  for  a  unit  or  combination  of  units  can  be  reduced  by  elim¬ 
inating  from  the  equation  (F)  =  (K)  •  (U)  joints  which  are  not  part  of  units  yet  to  be  considered. 
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The  new  matrix  will  represent  the  same  structure  as  the  original  matrix.  This  operation, 
which  amounts  to  the  partitioning  of  the  original  stiffness  matrix,  is  equivalent  to  the  elimina¬ 
tion  process  used  to  reduce  the  number  of  unknowns  in  a  set  of  simultaneous  equations  such  as 
Eq.  (2),  and  is  referred  to  as  "Elimination."  The  new  matrix  is  called  an  "Effective  Matrix.” 

To  illustrate  the  elimination  process  consider  the  force-displacement  equation  for  unit  01 
in  Fig.  2; 


“ 

• 

- 

*^0101 

'^OlOl 

^0102 

^0102 

*^0103 

^0103 

*^0104 

^0104 

^0105 

= 

K 

^0105 

*^0106 

^0106 

“^0304 

^0304 

*^0305 

^0305 

*^0306 

^0306 

where  F  and  U  represent  three  components  of  joint  force  and  joint  displacement,  respectively. 
Displacements  ^0103  eliminated  from  the  equations  (these  are  at  the 

joints  not  common  to  other  units)  and  Eq.  (4)  is  reduced  to 


“ 

ir< 

^0104 

^0104 

IT' 

^0105 

^0105 

*^0106 

K® 

^0106 

p» 

^^0304 

i\ 

^0304 

^0305 

'^0305 

pi 

^0306 

^0306 

The  original  matrix  (K)  of  size  27  x  27  is  thus  reduced  to  the  effective  matrix  (K^)  of  size  18  x  18. 
During  this  process,  the  force  vector  F  is  also  modified  to  F'.  The  details  of  this  modification 
are  described  in  Chapter  4. 

A  similar  effective  matrix  can  be  formed  for  unit  2  and  superimposed  on  that  for  unit  1.  The 
res\iltlng  matrix  represents  the  structure  shown  in  Fig.  5.  At  this  point,  the  fact  that  displace¬ 
ments  yQ|Q(  and  Sg^gg  must  equal  zero  is  taken  into  account  by  discarding  the  columns  and  rows 
in  the  matrix  corresponding  to  these  displacements.  The  matrix  for  the  two-imit  combination  is 
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then  25  X  25.  The  dieplacemente  at  joints  0105  and  0104  and  the  remaining  displacement  at  0106 
are  now  eliminated,  leaving  an  effective  16  x  18  matrix. 

Finally,  the  effective  matrix  for  unit  3  is  superimposed.  This  does  not  increase  the  size  of 
the  matrix  since  no  new  joints  are  introduced.  Joints  0204,  0205  and  0206  are  now  eliminated 
leaving  a  9  x  9  matrix.  The  displacements  at  joint  0306  are  all  zero  and  the  corresponding  terms 
in  the  matrix  may  be  discarded.  This  leaves  only  a  6  x  6  matrix,  which  is  finally  solved  for  the 
displacements  at  joints  0304  and  0305. 

By  the  simple  process  of  substitution,  all  other  displacements  may  be  determined.  Having 
displacements,  all  bar  stresses  and  reactions  are  easily  obtained. 

In  comparison  to  the  method  described  above,  a  direct  solution  without  subdivision  of  the 
structure  would  have  re,quired  the  solution  of  48  simultaneous  equations.  The  great  advantage 
of  the  method  is  now  obvious  when  applied  to  large  structures  having  many  joints. 


C.  INPUT  DATA 

It  is  desirable  to  prepare  the  input  data  on  format  input  data  sheets  before  punching  them 
onto  cards.  The  restrictions  of  the  STAIR  ANALYSIS  PROGRAM  and  the  items  of  data  required 
for  each  step  of  a  problem  are  listed  and  discussed  in  the  following  sections. 

Two  decks  of  cards  have  to  be  prepared  for  the  STAIR  System;  a  structural  data  deck  and 
a  load  data  deck.  These  decks  are  handled  differently  during  the  running  of  the  problem. 

The  structural  data  deck  is  input  directly  to  the  STAIR  System.  This  deck  begins  with  a 
block  of  cards  punched  from  an  "Initial  Data  Sheet."  Subsequent  blocks  of  cards  correspond  to 
consecutive  steps  in  the  ANALYSIS  PROGRAM;  the  first  card  of  each  block  calls  the  routine 
which  executes  the  step,  and  the  following  cards  provide  the  data  needed  for  the  execution. 

The  load  data  deck  is  input  to  the  STAIR  Load  Program  (SLOP),  This  program,  which  must 
be  run  before  the  ANALYSIS  PROGRAM,  receives  and/or  computes  concentrated  joint  loads  for 
various  types  of  loading  on  the  structure.  The  loads  are  assembled  into  groups  of  loads  applied 
to  individual  units  of  the  structure.  SLOP  is  stored  on  tape  for  use  by  the  STAIR  System. 

Each  item  of  data  entered  on  the  format  sheets  is  either  an  integer  number  (such  as  a  quan¬ 
tity  of  joints  or  bars)  or  a  fixed  number  (e.g.,  a  joint  coordinate  or  a  bar  area).  The  format 
sheets  prescribe  a  field  for  each  item  of  data.  A  field  is  the  maximum  niunber  of  characters 
an  item  may  contain. 

Integer  numbers  are  always  entered  at  the  right  of  their  field  and  are  never  written  with  a 
decimal  point.  Fixed  numbers  may  be  entered  in  one  of  two  ways:  with  or  without  a  decimal 
point.  If  the  decimal  point  is  specified  as  a  character  in  the  item  of  data,  the  number  may  be 
written  anywhere  within  its  field.  If  the  decimal  point  is  omitted,  the  number  must  be  positioned 
about  the  dashed  line  in  its  field,  with  the  fractional  part  to  the  right  of  the  dashed  line.  The  lo¬ 
cation  of  dashed  lines  is  written  into  the  machine  program  as  indicated  on  the  data  sheets  con¬ 
tained  herein.  If  there  is  no  dashed  line  in  the  field  of  a  fixed  number,  it  may  be  assumed  at 
the  right  of  the  field. 
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Minus  signs  must  always  be  specified  for  negative  numbers  and  must  always  be  written  im¬ 
mediately  before  the  first  nonzero  digit  in  the  field.  Plus  signs  need  not  be  written,  since  any 
unsigned  number  is  assumed  positive. 

The  following  examples  illustrate  these  points; 


Integer  Number  Format  Sheet  Read  by  Computer 


♦  I — I— 1— J — I  is  the  field  for  the  particular  items  whose  maximum  number  of  characters  is  four. 

D.  INPUT  DATA  FOR  THE  STAIR  LOAD  PROGRAM  (SLOP) 

The  STAIR  Load  Program  (SLOP)  is  written  independently  of  the  STAIR  System.  It  enables 
one  to  analyze  a  structure  subjected  to  more  than  one  loading  condition  without  repeating  the  com¬ 
plete  process.  The  input  data  of  SLOP  Is  stored  on  tape  and  is  ready  to  be  called  for  analysis 
by  the  STAIR  System.  The  data  sheets,  which  are  prepared  by  the  structural  analyst,  are  given 
on  pp.  14  to  29. 

Two  types  (INTP)t  of  loading  are  assembled  by  SLOP: 

Type  (1)  -  automatically  computed  joint  loads,  e.g.,  dead  load. 

Type  (2)  -  manual  input  joint  loads,  e.g.,  live  load. 

Each  loading  condition  may  be  Type  (1)  or  Type  (2)  or  a  combination  of  both  types.  Each 
condition  is  identified  by  an  NLD  (loading  condition  number);  if  the  same  NLD  number  is  given 
for  both  a  dead  and  a  live  load  condition,  the  analysis  will  be  made  for  the  combined  effect  of 
these  loads.  For  example,  on  pp.  19  and  20,  NLD  =  2  is  a  combined  dead  and  live  condition,  and 
NLD  =  3  or  4  are  conditions  of  live  load  only. 

Loading  condition  Type  (1)  may  be  computed  by  SLOP  for  any  orientation  of  the  structure; 
the  direction  of  the  pull  of  gravity  with  respect  to  the  coordinate  system  of  the  structure  is  spec¬ 
ified  by  Gx,  Gy  and  Gz,  the  direction  cosines  of  the  gravity  vector.  The  coordinates  of  the  cen¬ 
ter  of  gravity  and  the  total  weight  of  the  structure  are  output  from  each  dead  load  computation. 
The  computed  joint  dead  loads  are  stored  on  tape  for  the  STAIR  System. 

The  input  data  is  composed  of  four  groups  of  cards; 

1.  Initial  data  (refer  to  sheet  A,  p.l4). 

2.  Structural  data  (refer  to  sheet  B,  pp.  IS  to  17). 

3.  Loading  data  (refer  to  sheet  C,  pp.lB  to  22). 

(a)  Dead  load 

(b)  Live  load 


tSee  Loading  Data,  pp.18  to  22. 
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4.  Write-tape  data  (refer  to  sheet  D,  p.  23). 


BARSOL  data  and  some  of  the  STRIX  data  can  be  obtained  by  duplicating  the  input  data  of 
Group  2  in  the  SLOP  System.  SLOP  is  stored  on  tape  separately  from  the  STAIR  System,  and 
structural  data  are  required  for  dead  load  calculation  and  live  load  assembly;  therefore,  this 
input  data  must  be  repeated. 


Initial  Data  (sheet  A,  p.  14) 

(a)  Identification  problem  number  (IDP),  for  example,  200. 

(b)  Number  of  different  loading  conditions  to  which  the  structure  is 
subjected  (NK). 

(c)  Number  of  units  into  which  the  structure  is  divided  (NUT). 

(d)  Number  of  total  joints  in  the  structure  (NJT).* 

(e)  Number  of  total  bars  in  the  structure  (NBAR).* 

(f)  Number  of  joints  prefixed  with  each  unit  number  (lUADR).^  For 
example,  each  unit  of  the  example  structure  in  Fig.  2  has  9  joints, 
but  only  6  joints  are  prefixed  with  the  unit  number,  i.e.,  for  unit  01, 
joints  0101  to  0106. 

NUT 

^  (lUADR)j  =  Total  number  of  joints  in  structure 
i=l 

Structural  Data  (sheet  B,  pp.  15  to  17) 

(a)  Joint  card. 

(b)  Bar  card. 

Loading  Data 

Dead  Load  (sheet  C,  p.  18). 

(a)  INTP  =  1.  INTP  is  used  to  distinguish  the  following  three  blocks 
of  cards;  INTP  =  1  for  dead  load  input  block;  INTP  =  2  for  live 
load  block;  and  INTP  =  3  for  write-tape  input  block,  i.e.,  this 
instructs  ^e  computer  to  write  loads  on  tape  in  the  order  listed 
on  the  cards. 

(b)  Loading  condition  number  (NLD) . 

(c)  Density  of  material  used  for  the  structure  (DENS),  in  lbs.  per  cu.  ft. 

(d)  Direction  cosines  of  the  pull  of  gravity  relative  to  the  coordinate 
system  used  for  structure  (Gx,  Gy  and  Gz) . 

Note:  Items  (a),  (b),  (c)  and  (d)  are  repeated  for  each  loading  con¬ 
dition  containing  dead  load. 

Live  Load  (sheet  C,  pp.  20  to  22). 

(a)  INTP  =  2. 

(b)  Loading  condition  number  (NLD). 

(c)  Total  number  of  joints  loaded  (NJLD). 


♦  This  is  used  to  read  in  the  correct  niunber  of  joints  and  bars.  NJT  and  NBAR  can  be  zero  if 
no  dead  load  calculations  are  to  be  made. 

t Joints  are  always  consecutively  numbered  by  the  computer.  However,  if  the  programmer  num¬ 
bers  joints  0101,  0103,  0105,  0107,  then  lUADR  =  7,  not  4.  The  displacement  of  joints  0102, 
0104  and  0106  will  all  be  zero. 
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(d)  Joint  number  with  its  joint  load  in  components  (Jt.  No.,  Fx,  Fy 
and  Fz). 

Note:  Items  (a),  (b),  (c)  and  (d)  are  repeated  for  each  loading  con¬ 
dition  containing  live  load. 


Write -Tape  Data  (sheet  D,  p.  23) 


(a)  INTP  =  3. 

(b)  Unit  number  of  a  unit  (NUCK). 

(c)  Number  of  joints  in  the  unit  (NJU). 


(d)  Joint  card  giving  the  joint  numbers  contained  in  the  unit  in  row 
sequence.  Make  sure  the  joint  numbers  are  in  the  same  se¬ 
quence  as  the  joint  numbers  in  STRIX. 


Note:  Item  (a)  is  used  once,  whereas  items  (b) ,  (c)  and  (d)  are 
repeated  for  each  unit.  List  all  units  of  the  structure  in 
sequence. 


Size  Restrictions 

NK  =  number  of  loadings. 

NUT  =  number  of  units. 

NJT  =  total  number  of  joints  where  a  dead  load  is  to  be 
computed. 

NJU  =  number  of  joints  in  any  unit. 

NB  =  number  of  bars  input  In  this  step. 
lUAT  =  total  number  of  joints  in  structure. 

(a)  NUT  <  99. 

(b)  NB  <  1800.* 

(c)  NJT  <  800. t 

(d)  3  X  lUAT  X  NK  <  19000. 

(e)  3  X  NJU  X  NK<  2000. 


Output:  Dead  Load  Computation 

(a)  Coordinates  of  center  of  gravity. 

(b)  Total  weight  of  structure. 

Notes:  1.  These  outputs  should  be  the  same  for  any  orientation 
of  the  structure,  unless  the  density  is  changed. 

2.  The  joint  loads  are  assembled  on  tape  A8  in  binary 
form  for  input  to  STAIR. 

3.  If  Sense  Switch  3  on  the  computer  console  is  down, 
the  joint  loads  assembled  by  the  program  are  stored 
on  the  output  tape  (B2)  for  listing  (print  out). 

A  listing  of  the  input  deck  for  SLOP  for  this  problem  can  be  found  on  p.  25.  The  on-line 
print  out  and  the  output  of  the  SLOP  are  shown  in  pp.  27  to  29. 


'"NB  is  a  restriction  on  the  number  of  bars  whose  weight  can  be  calculated  in  any  one  SLOP  run. 
It  only  applies  if  gravity  loads  are  to  be  computed. 

tNJT  Is  a  restriction  on  the  number  of  joint  coordinates  that  can  be  stored  for  any  one  SLOP  run. 
It  only  applies  if  gravity  loads  are  to  be  computed. 
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28  LINES 


PROBLEM  2  on 


SLOP  INPUT  DATA  SHEET  *8' 
STRUCTURAL  DATA  II 


LIST  BAR  NUMBERS  OF  THE 

STRUCTURE  {2I5,F12.6) 

BAR  NUMBER 

AREA  (SQ.  IN) 

n 

IQ 

IB 

IB 

B 

IQ 

IB 

Id 

IH 

in 

If 

(E 

IIS 

m 

iE 

R 

in 

m 

\m\ 

L 

ir 

in 

in 

in 

in 

IB 

11 

in 

in 

a 

r 

■ 

m 

in 

in 

in 

IR 

IB 

in 

IP 

in 

a 

a 

■ 

in 

in 

in 

in 

in 

in 

in 

in 

in 

a 

r 

■ 

m 

n 

IB 

in 

■ 

in 

B 

ifi 

in 

in 

I 

■ 

i 

m 

III 

in 

in 

■ 

IR 

IB 

in 

in 

in 

a 

■ 

n 

n 

in 

in 

■ 

IR 

n 

in 

in 

n 

a 

■ 

m 

n 

in 

in 

in 

n 

IB 

n 

p 

a 

■ 

in 

n 

in 

in 

IR 

II 

IR 

in 

■ 

n 

a 

■nnn 

in 

■nnniniHi 

n 

a 

n 

■ 

n 

n 

in 

in 

■ 

R 

n 

R 

n 

:■ 

ET 

a 

■ 

m 

n 

n 

in 

iR 

n 

in 

in 

■ 

rn 

a 

■ 

t 

m 

n 

n 

in 

■ 

R 

n 

in 

n 

' 

a 

■ 

m 

13 

n 

in 

■ 

in 

R 

ifi 

n 

■ 

riS 

a 

n 

f] 

n 

in 

R 

n 

in 

na 

n 

n 

n 

in 

n 

n 

n 

n 

na 

n 

n 

n 

in 

n 

n 

n 

na 

■ 

n 

fi 

n 

in 

■ 

n 

n 

n 

■ 

■ 

na 

■ 

ra 

n 

n 

n 

■ 

n 

n 

n 

a 

■ 

■ 

na 

■ 

n 

n 

n 

n 

■ 

R 

n 

n 

n 

■ 

na 

r 

n 

n 

p 

p 

n 

n 

n 

n 

na 

■ 

n 

n 

n 

R 

■ 

R 

n 

R 

n 

■ 

■ 

m 

a 

■ 

n 

n 

n 

rs 

P 

■ 

r 

n 

R 

F 

■ 

■ 

na 

■ 

■ 

n 

|] 

n 

■ 

R 

n 

n 

n 

1 

na 

□ 

n 

n 

n 

n 

R 

n 

n 

n 

na 

■ 

n 

n 

n 

n 

■ 

n 

n 

n 

■ 

■ 

na 

a 

a 

r 

n 

n 

n 

r 

n 

R 

n 

na 

■ 

n 

n 

n 

n 

n 

n 

n 

n 

na 

r 

n 

n 

n 

RI 

R 

R 

n 

n 

rni 

n 

t 

n 

n 

R 

RI 

n 

n 

n 

n 

■ 

na 

a 

a 

n 

n 

n 

n 

R 

n 

n 

n 

r 

na 

a 

ri 

n 

n 

ni 

n 

R 

n 

R 

■ 

nn 

a 

j 

ni 

n 

ni 

R 

n 

R 

n 

■ 

1 

na 

a 

ni 

n 

R 

ni 

n 

n 

n 

ri 

2 

na 

:: 

! 

ni 

ni 

n 

ni 

p 

n 

n 

1 

! 

na 

n 

ri 

n 

ni 

r 

n 

n 

n 

na 

a 

ri 

n 

n 

Z]\ 

R 

n 

n 

n 

1 

mi 

PI 

n' 

Tl 

ni 

RI 

n 

1 

n 

! 

1 

nai 

pj 

ni 

n 

ni 

PI 

n' 

11 

n 

1 

mi 

ni 

ni 

R 

rii 

ni 

ni 

ni 

ni 

1 

1 

1 

mi 

□ 

ni 

ni 

□ 

rii 

■n 

ri 

nnm 

□ 

3 

□ 

□ 

mi 

im 

40UNES 


PAGE  I 

PUNCH  ONE  LINE 
PER  CARD- 
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PROBLEM  200 


SLOP  INPUT  DATA  SHEET  *8’ 
STRUCTURAL  DATA  II 


LIST  BAR  NUMBERS  OP  THE 

STRUCTURE  (21S,F)3.6) 

BAR  NUMBER 

AREA  (SO.  IN) 

n 

H 

IB 

IB 

B 

in 

fl 

Id 

m 

1 

m 

|[E 

'E 

IS 

IE 

m 

m 

a 

131 

1^1 

c_ 

ra 

1 

in 

in 

■ 

IR 

FI 

in 

n 

n 

■ 

II 

II 

■ 

in 

■ 

■ 

n 

n 

in 

n 

P 

ri 

m 

n 

II 

m 

II 

II 

in 

r 

n 

in 

in 

n 

n 

in 

in 

II 

■ 

II 

II 

in 

r 

R 

in 

in 

ir 

r 

n 

in 

r 

r 

r 

r 

n 

r 

r 

in 

n 

R 

n 

in 

in 

in 

r 

n 

n 

n 

■ 

r 

n 

IR 

IR 

n 

n 

r 

R 

n 

n 

n 

n 

IR 

n 

FT 

J 

■ 

R 

n 

n 

n 

■ 

p 

r 

in 

IG 

■ 

r 

R 

,n 

n 

n 

■ 

R 

ri 

in 

n 

n 

■ 

F 

n 

n 

R 

n 

n 

n 

n 

1 

FI 

n 

n 

n 

R 

m 

n 

n 

rn 

r 

n 

n 

R 

■ 

R 

n 

n 

n 

ni 

■ 

R 

n 

n 

R 

r 

n 

n 

n 

1 

m 

■ 

R 

R 

n 

■ 

R 

n 

■ 

ra 

1 

r 

■ 

■ 

n 

1 

L 

1 

n 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

n 

■ 

■ 

■ 

■ 

■ 

■ 

n 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

n 

c 

n 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

n 

l: 

_ 

n 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

n 

n 

■ 

■ 

■ 

a 

n 

■ 

■ 

■ 

■ 

■ 

n 

n 

L 

■ 

■ 

n 

■ 

■ 

■ 

■ 

■ 

■ 

n 

_ 

n 

■ 

■ 

■ 

■ 

■ 

n 

u. 

n 

n 

■ 

■ 

1 

n 

■ 

■ 

1 

1 

n 

L_ 

L_ 

■ 

■ 

■ 

■ 

■ 

■ 

1 

n 

u 

L 

□ 

□ 

□ 

□ 

■ 

■ 

I 

n 

□ 

J 

□ 

■ 

□ 

□ 

□ 

H 

: 

■ 

■ 

■ 

■ 

■ 

■1 

■ 

■1 

1 

■ 

n 

L 

□ 

□ 

□ 

□ 

■ 

■1 

■ 

■1 

1 

■ 

n 

J 

□ 

□ 

□ 

J 

■n 

L. 

Xi 
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PER  CARD 
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I 


SLOP  INPUT  DATA  SHEET  ‘C’ 


WH-  wifftif 


DATA  SHEET  'C 


PAGE  i/„ 


PUNCH  ONE  UNE  PER  CARO 


SLOP  INPUT  DATA  SHEET  ‘C* 
LOADING  DATA 
CASE'tP  (ShJ  efS"y 


20UNES 


20UNES 


SLOP  INPUT  DATA  SHEET  *€' 
LOADING  DATA 
CASEM  (Sh5'of^) 


20UNES 


PROBLEM 


SLOP  WRITE-TAPE  INPUT  DATA  SHEET  D 


PAGE  /4'n 


SLOP  WRITE-TAPE  INPUT  DATA  SHEET  D 

PROBLEM  Zoo 


PAGE 


INTP 

21 

MM 

begl 

oper 

•d  only  once  at 

nninflofthi,  PER  CARD 

ation  (II) 

(212) 

(715) 

_l _ 

Hi 

1  1 
NuckNJU 

_L 

4 

I. 

L 

1 

5. 

1£ 

11 

13 

14 

If 

Id 

17 

18 

19 

2( 

22 

22 

24 

27 

28 

29 

3C 

31 

32 

33 

34 

35 

r 

r 

r 

r 

r 

r 

r 

r 

r 

_ 

r 

n 

r 

2 

1-4 

(212) 

(715) 

Ha 

Nucl 

HI 

NJL 

1, 

2 

2} 

4, 

5 

71 

1 

9 

10 

n 

12 

K 

li 

15 

16 

17 

18 

19 

20 

21 

22 

24 

2S 

2d 

27 

28 

29 

3C 

32 

33 

33 

£i 

0. 

J_ 

Jl 

J 

St 

2. 

0 

a 

a 

i. 

i. 

o 

£ 

SL 

a 

4L 

k 

aj 

n 

£ 

1 

S. 

J! 

A 

6. 

k 

J 

J 

J 

•M 

mmm 

mmmm 

~1 

□ 

□ 

mtm 

■Ml 

. 

.... 

NOTE:  ORDER  OF  THE  JOINTS  MUST  BE  THE  SAME 
AS  IN  STRIX 
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200040300180054 

60606 


101 

00  . 

00 

0 

102 

00 

-3464102 

20 

103 

00 

-3464102 

-20 

104 

-60 

00 

00 

105 

-120 

-3464102 

-20 

106 

-120 

-3464102 

20 

201 

-30 

5196152 

00 

202 

-60 

6928203 

20 

203 

-60 

6928203 

-20 

204 

00 

10392305 

00 

205 

00 

17320508 

-20 

206 

00 

17320508 

20 

301 

30 

5196152 

00 

302 

60 

6928203 

20 

303 

60 

6928203 

-20 

304 

60 

00 

00 

305 

120 

-3464102 

-20 

306 

120 

-3464102 

20 

0101 

0102 

30 

0101 

0103 

30 

0101 

0104 

30 

0101 

0304 

30 

0102 

0103 

30 

0102 

0104 

30 

0102 

0105 

30 

0102 

0106 

30 

0102 

0304 

30 

0102 

0306 

30 

0103 

0104 

30 

0103 

0105 

30 

0103 

0304 

30 

0103 

0305 

30 

0103 

0306 

30 

0104 

0105 

30 

0104 

0106 

30 

0104 

0201 

30 

0104 

0202 

30 

0104 

0203 

30 

0105 

0106 

30 

0105 

0203 

30 

0106 

0203 

30 

0106 

0202 

30 

0201 

0202 

30 

0201 

0203 

30 

0201 

0204 

30 

0202 

0203 

30 

0202 

0204 

30 

0202 

0205 

30 

0202 

0206 

30 

0203 

0204 

30 

0203 

0205 

30 

0204 

0205 

30 

0204 

0206 

30 

0204 

0301 

30 

0204 

0302 

30 

0204 

0303 

30 

0205 

0206 

30 

0205 

0303 

30 

0206 

0302 

30 

0206 

0303 

30 

0301 

0302 

30 

0301 

0303 

30 

0301 

0304 

30 

0302 

0303 

30 

0302 

0304 

30 

0302 

0305 

30 

0302 

0306 

30 

0303 

0304 

30 

0303 

0305 

30 

0304 

0305 

30 

0304 

0306 

30 

0305 

0306 

30 

101 

1700  0 

102 

1700  0 

202 

0002 

202  5 

303  5 

203 
0005 

103  0 

202  5 

204  0 

206  0 

303  5 

204 
0001 

101  1 

3 


0 

0 


2 


-1 

-1 


75 

75 


0 

0 

0 

-75 

1 


3 


109 

0101  0102  0103  0104  0105  0106  0304 

0305  0306 

209 

0201  0202  0203  0204  0205  0206  0104 

0105  0106 

309 

0301  0302  0303  0304  0305  0306  0204 
0205  0206 


0 

0 


102 
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M604  STAIR-  SL0P2  LOAD  ASSEMBLY  FOR  PROBLEM  200 

FOR  INPUT  BLOCK  1*  TOTAL  WEIGHT  OF  STRUCTURE  IS  1»536  KIPS 

COMPONENTS  OF  THE  GRAVITY  VECTOR  ARE  0.  -1*00000  0* 

COORDINATES  OF  THE  CENTER  OF  GRAVITY  ARE  0.00000  3.46410  -O.OOOOO 

M504  STAIR-  SL0P2  LOAD  ASSEMBLY  FOR  PROBLEM  200 

FOR  INPUT  BLOCK  2*  TOTAL  WEIGHT  OF  STRUCTURE  IS  1*536  KIPS 

COMPONENTS  OF  THE  GRAVITY  VECTOR  ARE  0.  -1*00000  0* 

COORDINATES  OF  THE  CENTER  OF  GRAVITY  ARE  0*00000  3.46410  -0.00000 

M504  STAIR-  SLOP2  LOAD  ASSEMBLY  FOR  PROBLEM  200 

complete 


M504  SrfilR-  SLOF'2 

LOAD  ASSEMBLY 

FOR  PROBLEM 

Lonos 

FOR  UNIT  1 . 

LOADING  1 

101 

0. 

-0.035417 

0. 

102 

0. 

-0.104605 

0. 

103 

0. 

-0.104605 

0. 

104 

0. 

-0.097868 

0. 

105 

0. 

-0.084752 

0. 

10(S 

0. 

-0.084752 

0. 

304 

0. 

-0.097868 

0. 

305 

0. 

-0.084752 

0. 

306 

0. 

-0.084752 

0. 

LOflOS 

FOR  UNIT  1  , 

LOADING  2 

lot 

0. 

-0.035417 

0. 

102 

0. 

-0. 104605 

0. 

103 

0. 

-0. 104605 

0. 

104 

0. 

-0.097868 

0. 

105 

0. 

-0.084752 

0. 

106 

0. 

-0.084752 

0. 

304 

0. 

-0.097868 

0. 

305 

0. 

-0.084752 

0. 

306 

0. 

-0.084752 

0. 

LOfiOS 

FOR  UNIT  1, 

.  LOADING  3 

101 

0. 

0. 

0. 

102 

0. 

0. 

0. 

103 

0. 

-5.000000 

c . 

104 

0. 

0. 

0. 

105 

0. 

0 . 

0. 

106 

0. 

0. 

0. 

304 

0. 

0. 

0. 

305 

0. 

0. 

0. 

306 

0. 

0. 

0. 

LOf^tOS 

FOR  UNIT  1 , 

.  LOADING  4 

101 

1 .000000 

2.000000 

3.000000 

102 

0. 

0. 

0. 

103 

0. 

0. 

0. 

104 

0.. 

0. 

0. 

105 

0. 

0. 

0. 

106 

0. 

0. 

0. 

304 

0. 

0. 

0. 

305 

0. 

0. 

0. 

306 

0. 

0. 

0. 

LOfiOS 

FOR  UNIT  2 

.  LOADING  1 

201 

0. 

-0.035417 

0. 

202 

0. 

-0. 1 04 605 

0 . 

203 

0. 

-0 . 1 04605 

0. 

204 

0. 

-0.097868 

0. 

205 

0. 

-0 . 084752 

0. 

206 

0. 

-0 . 084752 

0  , 

104 

0. 

0. 

0 . 

105 

0. 

0. 

0 . 

106 

0. 

0. 

0. 

LOfiOS 

FOR  UNIT  2 

.  LOADING 

201 

0. 

-0.035417 

0. 

202 

5.000000 

-0.1 04605 

7.500000 

203 

0. 

-0. 1 04605 

0. 

204 

0. 

-0 . 097868 

0. 

205 

0. 

-0.084752 

0. 

206 

0. 

-0.084752 

0. 

104 

0. 

0. 

0. 

105 

0. 

0. 

0. 

106 

0. 

0. 

0. 

200 
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r 
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LOfiOS  FOR  UMIT  2. 

loading  3 

201 

0. 

0. 

202 

5.000000 

0. 

203 

0. 

0. 

204 

0. 

•10.000000 

205 

0. 

0. 

206 

0. 

0. 

lO'* 

0. 

0. 

105 

0. 

0. 

106 

0. 

0. 

LOADS 

FOR  UNIT  2. 

LORDING  4 

201 

0. 

0. 

202 

0. 

0. 

203 

0. 

0. 

204 

0. 

0. 

205 

0. 

0. 

206 

0. 

0 . 

104 

0. 

0. 

105 

0. 

0. 

106 

0. 

0. 

LOADS 

FOR  UNIT  3. 

LOADING  1 

.301 

0. 

-0.035417 

302 

0. 

-0.104605 

303 

0. 

-0.104605 

304 

0. 

0. 

305 

0. 

0. 

306 

0. 

0 . 

204 

0. 

0. 

205 

0. 

0 . 

206 

0. 

0. 

LOADS 

FOR  UNIT  3. 

loading  2 

301 

0. 

-0.03541 7 

302 

0. 

-0. 1 04605 

303 

5.000000 

-0. 104605 

304 

0. 

0 . 

305 

0. 

0. 

306 

0. 

0. 

204 

0. 

0. 

205 

0. 

0 . 

206 

0. 

0. 

LOADS 

FOR  UNIT  3; 

,  LOADING  3 

301 

0. 

C- 

0 . 

0  . 

302 

303 

•J  • 

5. CO 0000 

0. 

304 

0. 

0. 

305 

0. 

0. 

306 

0. 

0. 

204 

0. 

0. 

205 

0. 

0. 

206 

0. 

0. 

LOADS  FOR  UNIT  3 

,  LOADING  4 

301 

0. 

0. 

302 

0. 

0. 

303 

0. 

0. 

304 

0. 

0 . 

305 

0. 

0. 

306 

0. 

0. 

204 

0. 

0. 

205 

0. 

0. 

206 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

-7.500000 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

-7.500000 

0 . 

0. 

0. 

0. 

0. 

0. 


0 . 

0. 

1 .000000 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


E.  THE  STAIR  OPERATIONS 


The  STAIR  System  is  divided  into  ten  different  operations; 
Code  No.  and  Name*  Operation 


1 


2 


3 


1 


5 


I 


10 


STRIX 


EFFRIX 


MATRO 


MADD 


MATOUT 


MATIN 


DEFSOL 


BARSOL 


CHECK 


STOP 


Forms  a  x^iffness  matrix  from  an  input  of  joint 
coordinates,  bar  areas  and  elastic  modulus. 

Eliminates  a  joint  or  a  group  of  joints  from  a 
stt^ness  mat^  equation  and  the  row  sequence 
which  Identifies  the  stiffness  matrix. 

Transforms  the  coordinates  of  a  stiffness  matrix  to 
give  a  new  matrix  of  another  identical  unit.  This 
transformation  consists  of  the  station  of  two  axes 
about  the  third  axis. 

(Matrix  addition) 

Superimposes  two  stiffness  matrices:  one  in  the  computer 
core  and  the  other  on  tape. 

(Matrix  out) 

Transfers  a  matrix  from  the  operating  core  of  the 
computer  to  magnetic  tape. 

(Matrix  in) 

Transfers  a  matrix  stored  on  magnetic  tape  back  into 
the  core  of  the  computer. 

(Deflection  solution) 

Computes  joint  displacements  from  a  stiffness  matrix 
which  completely  represents  the  structure  under  a 
given  set  of  loading  conditions. 

(Bu*  stress  solution) 

Computes  the  stress  in  any  bar  from  the  displacements 
of  its  end. 

Sums  all  the  bar  forces  at  each  joint  in  each  coordinate 
direction. 

Ends  the  analysis  sequence. 


'''Whenever  an  operation  is  given  on  a  data  sheet,  t>oth  the  code  number  and  name 
must  appear. 

Any  three-dimensional  pin-jointed  structure  can  be  analyzed  by  certain  combinations  of  these 
operations.  Such  an  operation  sequence  is  called  an  ANALYSIS  PROGRAM.  The  particular 
ANALYSIS  PROGRAM  used  should  perform  the  required  analysis  in  a  minimum  of  computer  time 
without  exceeding  the  size  limitation  on  matrices  imposed  by  the  computer  capacity.  The  com¬ 
puter  time  required  for  each  operation  is  different.  For  the  details  of  computer  time  accumula¬ 
tion,  an  experienced  progranuner  should  be  consulted. 


The  elze  of  a  stlffnesa  matrix  in  operation  is  limited  by  the  capacity  of  the  computer  uaed. 
The  size  of  matrix  which  represents  a  structure  increases  with  the  square  of  three  times  the 
number  of  joints;  the  largest  number  of  joints  which  can  be  handled  by  an  ANALYSIS  PROGRAM 
at  any  time  in  the  IBM-7090  computer  is  59.  This  is  the  main  reason  for  dividing  the  structure 
into  units  and  is  also  why  EFFRDC  is  used  to  eliminate  the  joints  (or  reduce  the  size  of  matrix) 
as  soon  as  possible  in  the  ANALYSIS  PROGRAM. 

Several  other  parameters  are  limited  in  magnitude  in  the  STAIR  System.  These  limitations 
are  listed  in  the  next  section.  The  STAIR  System  detects  violations  of  these  restrictions  and 
prints  an  Indication  of  the  error. 

A  STAIR  ANALYSIS  PROGRAM  SHEET  with  a  routine  to  analyze  the  example  structure  in 
Fig.  2  is  presented.  This  sheet  should  be  completed  before  the  formation  of  input  data  for  each 
step  of  the  analysis.  The  following  symbols  are  used  in  this  description: 


[  ]  contained  in  computer  core. 

(  )  contained  in  tape  location  -  1. 
((  ))  contained  in  tape  location  —  2. 
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Step  I  —  STRIX  ferim  lh«  stlffnM  matrix  for  unit  01  which  It  Idontlflod 
by  row  Mquonca  0101  —  0106,  0304  —  0906. 


Stop  2-  EFFRIX  •IlminotM  |olntt  0101  -  0103  from  unit  1. 


Stop  3  -  MATOUT  trantfara  th*  motrix  In  cor*  to  mognatle  top*  loc.  1. 


Stop  4-  MATRO  tranifonra  tha  matrix  In  cor*  Into  matrix  for  unit  2. 


St*p  5-  MADO  fup*rimpoMi  Ih*  matrix  in  cor*  and  matrix  in  tap*  loc.  I. 


Step  7  -  MATOUT  traiwfcrt  th«  matrix  In  eora  to  tap*  location  2  • 


c 


n 


•ffoetlva 
unit  1) 


Stop  8-  MATIN  tranttora  the  matrix  In  top#  locotlon  I  to  th«  eora. 


•ffoctiva 
unit  1) 


Stop  9-  MATRO  trantformi  the  matrix  in  eora  Into  matrix  for  unit  3. 


•ffoctiva 
unit  1) 


•ff.  unit  1 
•ff.  unit  2^^ 


( 


Slap  10  -  MAOO  mpfliwpem  tli*  matrix  In  eera  and  lha  matrix  In  topa 
location  2  In  tha  raw  laquanea  0204  -  0206/  0304  -  0306. 


afhctiva 
unit  1) 


Slap  12  -  OEFSOL  eomputat  the  joint  dltplacamanti  of  all  joinh  and  raoetlon*. 
Stop  13  -  BARSOL  computM  tha  bar  itrowat. 

Stop  14  -  CHECK  givar  lummatlon  of  Intamal  fbreat  at  aoeh  joint. 

Stop  15  -  STOP  ondi  tfia  anal)rth  program. 


F.  INPUT  DATA  FOR  THE  STAIR  SYSTEM 


The  Input  parameters  required  by  each  step  of  the  STAIR  operations  are  described  below. 
The  restrictions  on  the  magnitude  of  these  parameters  and  operations  are  given  where  appli¬ 
cable.  The  formats  In  which  these  data  parameters  are  prepared  will  be  illustrated  In  the  fol¬ 
lowing  data  sheets.  Data  sheets  for  each  step  are  arranged  In  the  same  order  as  the  example 
for  the  ANALYSIS  PROGRAM.  The  same  example  structure  of  Fig.  2  Is  used. 

Initial  Input  (refer  to  data  sheet,  p.  44) 

(a)  Problem  identification  number  (IDP) . 

(b)  Number  of  units  Into  which  the  given  structure  is  divided  (NUT). 

(c)  Number  of  different  loading  conditions  to  which  this  given  structure 
Is  subjected  (NK). 

(d)  Number  of  reaction  joints  or  supports  in  structure. 

(e)  Elastic  modulus  of  bar  material  E  kips  per  sq.  in. 

(f )  List  of  reaction  joint  numbers  and  the  coordinate  directions  in 
which  displacement  Is  restricted  (1-immovable;  O-free). 

Size  Restrictions 

(a)  Number  of  units  <  99. 

(b)  Niunber  of  reaction  joints  <  125. 

(c)  Loading  conditions  <  28. 

STRIX  (refer  to  data  sheet,  p.  45) 

(a)  Step  munber  in  ANALYSIS  PROGRAM.  It  precedes  every  operation 
code,  as  the  001  on  p.45  precedes  01  STRIX,  and  the  005  precedes 
04  MADD  on  p.  49.  It  is  the  first  card  of  a  card  block  calling  the 
routine  which  executes  the  step  as  mentioned  in  the  section  on  input 
data.  (Step  numbers  must  start  at  01  and  proceed  sequentially  to 
the  end.) 

(b)  Identification  number  of  th«*  unit  (NUU)  to  which  this  STRIX  operation 
applies.  On  p.  45  NUU  is  01;  later,  units  02  and  03  will  be  formed 
by  MATRO.  In  case  some  of  the  units  are  not  identical,  NUU  is  used 
to  identify  the  different  STRIX  operations  for  different  units  when 
STRIX  is  called  more  than  once. 

(c)  Number  of  joints  in  the  unit  (NJTU). 

(d)  Number  of  bars  in  the  unit  (NMU). 

(e)  Joint  sequence.  This  serves  as  a  reminder  only,  and  is  not  a  part 
of  the  input  data. 

(f)  Joint  cards  giving  the  joint  numbers  [see  Operation  Restriction  (b)] 
contained  in  the  unit  in  sequence*  and  the  coordinates  of  the  joints 
in  feet  (reproduced  from  SLOP). 

(g)  Bar  cards  are  reproduced  from  SLOP.  Bars  which  are  common  to 
two  units  are  listed  separately  for  this  input  (see  p.  46). 


*The  sequence  can  be  of  any  convenient  form,  but  the  joints  that  are  to  be  eliminated  by  EFFRIX 
must  be  placed  first  in  the  sequence. 


Size  Reatrictlons 


If  N  =  3x  (numbei*  of  joints  In  the  unit),  NK  =  number  of  loading  conditions, 

(a)  +  N  X  NK  <  16000  (see  table). 


N 

NK 

59 

1 

56 

4 

57 

7 

56 

10 

55 

13 

54 

17 

53 

20 

52 

24 

51 

27 

(b)  Number  of  bars  <  600. 

(c)  Unit  number  <  total  number  of  units  in  the  structure. 

Operation  Restrictions 

(a)  In  no  case  may  all  the  bars  meeting  at  a  joint  lie  in  one  plane  unless 
the  joint  la  a  reaction  joint.  If  a  planar  joint  exists,  a  new  bar  out 
of  the  plane  (an  Imaginary  member)  may  be  inserted  which  will  carry 
no  stress  but  will  allow  the  analysis  to  proceed.*  The  new  bar  may 
connect  the  planar  joint  with  any  other  joint  not  lying  in  the  plane. 

(b)  If  an  EFFRIX  operation  immediately  follows  STRIX  in  the  sequence 
of  the  ANALYSIS  PROGRAM,  those  joints  which  are  going  to  be  elim¬ 
inated  have  to  be  placed  first  ^  the  row  of  joint  cards.  In  the  ex- 
ample,  joints  0101,  6l62  and  0103  are  eliminated  by  the  following 
EFFRIX. 

EFFRIX  (refer  to  p.  47) 

(a)  Step  nximber  in  ANALYSIS  PROGRAM. 

(b)  Number  of  joints  (NE)  to  be  eliminated  from  the  row  sequence 
established  by  the  previous  operation  (STRIX  or  MADD) . 


Size  Restrictions 

NE  =  number  of  joints  to  be  eliminated. 

NJ  =  number  of  joints  in  original  row  sequence  (number 
of  joints  in  the  unit) . 

NR  =  number  of  reaction  joints  among  joints  being  eliminated. 

NC  =  munber  of  reaction  components  acting  at  the  eliminated 
reaction  joints. 


*  See  p.  41  for  a  more  detailed  explanation  of  a  planar  joint. 


NK  =  number  of  different  loading  conditions. 

A  =  3NE  -  NC. 

B  =  3(NJ~NE). 

C  =  3NE. 

(a)  NR<  to. 

(b)  +  A(B  +  NK)  <  9500. 

(c)  A^<(Cx  +  CXB1. 

(d)  <  [C  X  +  C  X  B]. 

MATOUT  (refer  to  p.  47) 

(a)  Step  number. 

(b)  Tape  location  (Loc.)  in  which  a  matrix  is  to  be  stored.  Matrices 
are  stored  on  tape  in  blocks.  These  are  numbered  consecutively, 
the  tape  location  being  the  number  given  to  the  block. 


Operation  Restriction 

Only  one  matrix  may  be  stored  in  a  tape  location. 


MATRO  (refer  to  p.  48) 


(a)  Step  number. 

(b)  Identification  nvunber  of  new  unit  (NUN). 

(c)  Angle  through  which  old  unit*  must  be  rotated  in  the  x-y  plane  to 
arrive  at  the  position  of  the  new  unit.  Angle  (EPS)  measured  in 
degrees  and  positive  when  clockwise,  as  viewed  by  an  observer 
on  the  positive  z-axis.  In  the  STAIR  system,  the  rotation  must 
be  about  the  z-axis. 

(d)  Row  sequence  of  the  resulting  matrix  after  rotation  [NROW  (1) 
to  NROW  (NJT)). 

(e)  If  the  matrix  being  rotated  is  a  reduced  matrix,  i.e.,  some  joints 
have  been  eliminated  from  it  by  EFFRIX  in  a  previous  step,  the 
inverse  matrix  corresponding  to  these  eliminated  joints  is  stored 
in  the  secondary  matrix  storage  -  then,  not  otherwise.  The  new 
row  sequence  of  the  eliminated  joints,  i.e.,  the  new  row  sequence 
for  the  new  inverse  matrix,  is  NROSEC  (1)  to  NROSEC  (NE). 

Size  Restriction 

(a)  New  unit  number  <  number  of  units  in  the  structure. 

Operation  Restriction 

(a)  Matrices  of  identical  units  from  which  joints  have  been  eliminated 
may  be  formed  by  MATRO  except  when  one  of  the  eliminated  joints 
of  the  old  unit  has  a  corresponding  joint  in  the  new  unit  which  is  a 
reaction  joint,  or  vice-versa. 

(b)  MATRO  can  not  follow  a  double  EFFRIX  sequence. 


♦Matrix  of  unit  in  core. 


MADD  (refer  to  p.49) 


(a)  Step  number. 

(b)  Number  of  joints  In  the  resulting  structure  (NF)  not  including  those 
previously  eliminated. 

(c)  Tape  location  (Loc.)  in  which  a  matrix  has  been  stored  by  MATOUT. 

The  matrix  in  this  tape  location  is  superimposed  on  the  matrix  in 
the  computer  core  by  this  operation. 

(d)  List  of  row  sequence  of  the  resulting  matrix.  Caution  should  be  taken 
in  setting  up  the  row  sequence  if  EFFRIX  is  planned  to  follow  the  MADD 
operation.  Those  joints  which  are  being  eliminated  in  EFFRIX  should 
be  placed  first  in  the  row  sequence. 

Size  Restriction 

NF  =  ntunber  of  joints  in  resulting  structure. 

NA  =  larger  number  of  joints  in  either  unit  being  added. 

NK  =  number  of  loading  conditions. 

(a)  3NF(3NF  +  l)/2  +  3NF  X  NK  <  16000  (see  table). 

(b)  3NA(3NA  +  l)/2  +  3NA  X  NK  <  9500  (see  table). 

(c)  NA  <  45. 


NF 

NK 

59 

1 

58 

4 

57 

7 

56 

10 

55 

13 

54 

17 

53 

20 

52 

24 

51 

27 

NA 

NK 

45 

2 

44 

5 

43 

8 

42 

11 

41 

15 

40 

18 

39 

22 

38 

25 

37 

28 

EFFRIX  (refer  to  p.  50,  similar  to  Step  2) 

MATOUT  (refer  to  p.  50,  similar  to  Step  3) 

MATIN  (refer  to  p.  50) 

(a)  Step  number. 

(b)  Tape  location  (Loc.)  of  the  matrix  to  be  read  into  computer  core. 

Operation  Restriction 

None. 


MATRO  (refer  to  p.  51,  similar  to  Step  4) 
MADD  (refer  to  p.  52,  similar  to  Step  5) 
EFFRIX  (refer  to  p.  53,  similar  to  Step  2) 
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DEFSOL  (refer  to  p.  54) 


t 


(a)  Step  number. 

(b)  List  of  the  number  of  joints  which  are  prefixed  by  each  unit  number 
(lUADR).  For  example,  there  are  9  joints  In  unit  01,  but  only  6 
joints  are  prefixed  by  unit  number  01.  (Notes  regarding  this  In 
SLOP  apply  here  also.) 


{ 


i 


Size  Restrictions 

NJ  =  number  of  joints  In  final  row  sequence. 

NR  =  number  of  reaction  joints  among  final  joints. 

NC  =  number  of  reaction  components  among  final  joints. 
NK  =  number  of  loading  conditions. 

NJT  =  total  number  of  joints  in  structure. 

NRT  =  total  number  of  reaction  joints  in  structure. 

A  =  3NJ  -  NC. 

(a)  NR  <  10. 

(b)  A  X  (A  +  l)/2  +  A  X  NK  <  9408. 

(c)  A^  <  16,000. 

(d)  4  .  NK(NJT  +  NRT)  <  16,000. 


Operation  Restriction 


DEFSOL  may  be  called  only  when  the  final  matrix  of  the  structure  is  stored  in  the  computer 


core. 


Output 


(a)  Joint  displacement  (in  inches). 

(b)  Reactions  (In  kips). 


BARSOL  (refer  to  p.  54) 


(a)  Step  number. 

(b)  Number  of  bars  whose  stresses  are  to  be  computed  (NM). 

(c)  Number  of  joints  to  which  the  bars  listed  in  (b)  are  connected  (NJ). 

(d)  List  of  bar  numbers  of  the  bars  whose  stresses  are  to  be  computed 
(reproduced  from  SIX)P). 

Note:  STAIR  will  compute  any  desired  number  of  bar  stresses  by 
repeated  use  of  BARSOL. 

(e)  Joint  cards,  giving  the  joint  numbers  and  coordinates.  Only  those 
joints  which  appear  in  the  list  of  bar  numbers  (p.  54)  need  be  in¬ 
cluded  (reproduced  from  SLOP) . 

Size  Restrictions 


NK  =  number  of  loading  conditions. 

(a)  NJ<  800. 

(b)  NM  <  1000. 

(c)  NKXNM<4300. 
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Operation  Reatrictlon 

BARSOL  may  be  called  only  after  the  DEFSOL  operation. 


Output 

(a)  Bar  stresses  (kips  per  sq.  in.)  for  all  loading  conditions. 

CHECK  (refer  to  p.55) 

(a)  Step  number. 

(b)  Total  number  of  bars  in  structure. 

(c)  Total  number  of  joints  in  structure. 


Size  Restriction 

None. 

Output 

Sum  of  all  b»  forces  at  each  joint  (in  kips)  for  all  loading  conditions. 

STOP  (refer  to  p.  55) 

(a)  Step  number. 


Notes  on  Input  Data  Units 

The  above  description  assumed  that  the  input  data  are  as  follows; 

Joint  coordinates  (x,  y,  z)  (in  ft) 

Bar  areas  (A)  (in  sq.  in.) 

Modulus  of  elasticity  (E)  (in  kips  per  sq.  in.) 

Loads  (F)  (in  kips) 

The  results  are  given  as  follows; 

Deflections  (in  in.) 

Reactions  (in  kips) 

Bar  stresses  (in  kips  per  sq.  in.) 

Equilibrium  check  (in  kips) 

The  input  data  could  be  in  any  units,  but  in  such  a  case,  the  output  would  be  in  terms  of  the 
following  units; 

FL 

The  deflections  would  be  in  the  units  of  ,  and  bar  stresses  would 
be  in  the  units  of  E, 

Reactions  would  be  in  terms  of  F,  and  equilibrium  check  in  terms  of  F. 

A  listing  of  the  input  cards  for  the  ANALYSIS  PROGRAM  is  shown  on  p.  56,  and  the  output  from 
a  successful  run  for  the  sample  problem  is  shown  on  p.59. 
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G.  STAIR  SYSTEM  ERROR  STOPS 


The  STAIR  system  detects  several  errors,  in  addition  to  violation  of  size  restrictions.  When 
an  error  is  detected,  the  problem  number  and  the  step  of  the  ANALYSIS  PROGRAM  at  which  the 
error  occurs  are  printed  by  the  computer,  as  well  as  a  diagnosis  of  the  trouble. 

Some  error  stops  may  be  bypassed  by  setting  Sense  Switch  1  down;  this  may  be  necessary 
under  unusual  circumstances,  but  care  should  be  exercised.  Depressing  Sense  Switch  1  on  the 
machine  console  causes  the  following  error  stops  to  be  bypassed; 

STRIX  error  (a)  (see  below), 

SLOP  errors  (a),  (b)  and  (d)  (see  below). 

In  each  case,  the  error  diagnosis  is  printed  before  the  program  continues. 

Some  errors  are  not  caused  by  mistakes  in  input,  but  are  machine  errors.  A  common  com¬ 
plaint  arises  from  the  malfunctioning  of  the  tape  units.  When  machine  error  is  suspected,  a 
qualified  programmer  should  be  consulted  before  rerunning  the  problem. 

Error  diagnoses  written  by  the  various  STAIR  routines  are  given  below. 

MISLAM* 

A  FORTRAN  II  MISLAM  diagnosis  indicates  a  machine  error  in  reading  the  STAIR  system 
tape.  It  may  be  necessary  to  rewrite  this  tape  from  the  binary  card  deck  containing  all  the 
STAIR  routines. 

CONTROL 

(a)  Illegal  subroutine  code  used;  incorrect  data  for  a  preceding  step  or 
a  mispunched  call  card. 

(b)  Input  data  file  out  of  order;  a  wrong  step  number  on  a  call  card  or 
incorrect  data  for  a  preceding  step. 

(c)  Size  limitations  are  exceeded;  mistake  on  data  sheet. 

(d)  Program  pause;  used  for  program  testing,  probably  indicates  an 
illegal  subroutine  code  (see  (a)  above]. 

STRIX 

(a)  The  following  joints  are  planar  joints;  either  the  wrong  coordinate  data 
are  supplied  or  the  bars  at  the  joints  listed  lie  in  a  plane  .t  If  there  is  no 
load  acting  normal  to  the  plane  of  the  bars,  a  new  bar  may  be  inserted 
which  will  carry  no  stress  but  which  will  allow  the  analysis  to  proceed. 

(b)  Space  limitations  have  been  exceeded;  check  input  parameters. 

(c)  Joint  N  used  in  a  bar  number,  but  is  not  in  the  row  sequence;  check 
bar  number  and  row  sequence. 

(d)  Loads  can  not  be  found;  tape  error.  (Tape  AS  should  contain  the  loads.) 


’’'See  Chapter  4  for  description  of  MISLAM. 

tThe  check  for  a  planar  joint  is  performed  in  the  following  manner;  Let  =  unit  vector  in  di¬ 
rection  of  bar  1;  L2  =  unit  vector  in  direction  of  bar  2;  L  =  unit  vector  in  direction  of  any 
other  bar.  First,  the  vector  product  C  =  Lj  k  L2  is  executed.  If  each  component  of  C  is  less 
than  0.000001,  the  bars  are  considered  collinear,  and  the  unit  vector  associated  withthe  next  bar 
is  used  as  L2  to  form  a  new  'C.  Then,  the  scalar  product  /C-L/  is  executed.  If  /C'L/  < 
0.00001  for  every  remaining  bar,  the  joint  is  considered  planar. 


effrix 


(a)  Space  limitations  have  been  exceeded;  check  input  parameters. 

(b)  Computer  stops  without  diagnosis;  possibly  the  matrix  is  singular 
(unstable  structure).  Check  the  previous  STRIX  stops  for  planar  joints 
now  being  EFFRIXed. 

MATRO 

(a)  An  attempt  has  been  made  to  rotate  the  effective  matrix  of  a  unit  with 
internal  reaction  joints;  the  ANALYSIS  PROGRAM  must  be  changed. 

Either  the  old  or  the  new  unit  contains  reactions  at  the  joints  which 
were  eliminated. 

(b)  The  number  of  the  new  unit  is  wrong;  check  input  to  this  step  and  the 
total  number  of  units  as  given  in  the  initial  input  block. 

(c)  Loads  can  not  be  found;  tape  error  (see  STRIX). 

MADD 

(a)  Space  limitations  have  been  exceeded;  check  input  parameters. 

(b)  Joint  N  is  not  included  in  the  final  row  sequence;  either  the  final  row 
sequence  input  at  this  step  is  wrong,  or  the  row  sequence  of  one  of  the 
matrices  being  added  is  wrong. 

(c)  Matrix  on  tape  can  not  be  found;  tape  error  (check  location  number). 

MATIN  or  MATOUT 

(a)  Record  group  can  not  be  found;  tape  error. 

DEFSOL 

(a)  Space  limitations  are  exceeded;  check  input. 

(b)  No  diagnosis;  see  EFFRIX  and  previous  comments  on  machine  errors. 

EARSQI, 

(a)  Space  limitations  have  been  exceeded;  check  input. 

(b)  Joint  coordinate  not  defined  for  this  bar,  the  coordinates  of  one  end 
of  the  bar  were  not  supplied,  the  bar  number  is  illegal  or  a  joint 
number  in  the  input  or  coordinates  is  illegal.  (This  error  is  printed 
in  the  output,  but  does  not  stop  the  program.) 

SLOP 

(a)  Size  limitations  are  exceeded;  check  input. 

(b)  Illegal  joint  numbers  or  bar  numbers  (list);  error  in  the  input  of  joint 
coordinates,  bar  numbers,  loads  or  unit  row  sequences. 

(c)  Illegal  input  heading;  an  input-type  code  has  been  read  which  is  not  1, 

2  or  3. 

(d)  Joint  list  out  of  consecutive  order;  row  sequence  of  units  is  out  of  order 
or  joint  number  is  wrong. 

Illegal  joint  numbers  in  BARSOL  and  SLOP  are  detected  when  the  first  two  digits  of  a  joint 
number  are  zero  or  greater  than  the  number  of  units,  or  when  the  last  two  digits  are  greater 
than  the  number  of  joints  in  the  unit  which  contains  the  joint. 
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PROBLEM  ?O0 


STAIR  INITIAL  DATA  SHEET 


PAGE  Z//3 


PUNCH  ONE  LINE  PER  CARD 
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PROBLEM  20  O 


PAGE  -4//3 


STAIR  STRIX  DATA  SHEET  'B' 


Bar  numbers  of  operation  STRIX  (and  of  any  following  STRIX  operation)  that  are  common 
to  adjacent  units  mutt  have  new  bar  cards  punched  using  1/2  the  actual  bar  area.  All 
other  bar  numbers  of  a  STRIX  operation  can  use  reproduced  bar  cards  taken  from  the 
SLOP  deck. 


(2I5>  FI 2.6) 


PUNCH  ONE 
LINE  PER  CARD 


The  order  of 
STRIX  bar  cards 
is  immaterial 
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STAIR  EFFRIX  DATA  SHEET 


PAGE  r//» 


PUNCH  ONE  LINE  PER  CARD 


STAIR  MATIN  DATA  SHEET 


PROBLEM 


PAGE  / 

PUNCH  ONE  LINE  PER  CARD 


3  5  _ y 

!  !'  ^^a!t  ■'  liNt  03,12, 6H) 
Step  N(  Operation 


(12) 


PROBLEM  Z60 


STAIR  MATRO  DATA  SHEET 


PAGE 
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STAIK  MAOD  DATA  SHEET 
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STAIR  MATIN  DATA  SHEET 


PROBLEM  ^00 


PAGE  / 

PUNCH  ONE  LINE  PER  CARD 


3  5  11 

o!o!8  m mIa!t  |!n!~  03,12, 6H) 

Step  Nc  Operation 


PROBLEM  loo 


STAIR  MATRO  DATA  SHEET 


PAGE  ^//3 


STAIR  MADD  DATA  SHEET 


PROBLEM 


PAGE 
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STAIR  EFFRIX  DATA  SHEET 

PROBLEM 0  PAGE  ////f 


53 


STAIR  DEFSOL  DATA  SHEET 


PROBLEM  O 


PAGE  / V/3 


3l  S|  111 


BcaatiajMijfco 


St«^No.  I  I  Op«p«»l«n  I 


a3,I2,iH) 


PUNCH  ONE  LINE  PER  CARO 


Ll«l  nuinb«r  el  jeintt  which  ore  prefixed  with  a 
unit.  Lilt  them  in  sequence  of  units. 

Punch  row-wise,  one  line  per  cord. 


PROBLEM 


STAIR  BARSOL  DATA  SHEET 


PAGE  / 


3  I  SI  H 

0»  B  iaIrIs  [oil 

Step  No.  Operetlen 


(I3,I2.6H) 


_ 

_ 101 

man 

(215) 

No.«f  B«r« 

No.  of  Jto. 

PUNCH  ONE  LINE  PER  CARD 


BARSOL  |oint  and  bar  cords  to  bo  reproduced  from  the 
SLOP cords.  The  order  of  cords  is  imnioteriol,  except 
area  cords  must  come  before  the  geometry  cords. 


PROBLEM  Zoo 


STAIR  CHECK  DATA  SHEET 


PAGE  /3//3 


STAIR  STOP  DATA  SHEET 


PROBLEM  ZOO  PAGE  / 
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1  IStRlX 

2  2EFFK1X 
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If  4FADO 
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12  TOEFSCL 

MSi4»  STAIR  PROBLEK  NUMBER  2»t 
DISPLACEMENTS  FOR  LOADING  CONDITION  1  FOLLOW 
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DISPLACEMENTS  FOR  LOADING 

CONDITION  2  FOLLOW 
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OISPLACEPENTS  FOR  LOADING 

CONDITION  3  FOLLOW 
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-0.00529 

V« 

0.01051 

w* 

0.12555 

JOINT 

203 

u* 

0.00221 

V- 

0.04171 

w* 

0.12564 

JOINT 

204 

u« 

0.00604 

V* 

0.01570 

w* 

0.15656 

JOINT 

205 

u* 

0.01374 

V« 

0.03230 

w* 

0.20789 

JOINT 

206 

u« 

0.01927 

V« 

-0.00131 

w« 

0.20753 

JOINT 

301 

u» 

0.00401 

V» 

0.01655 

w« 

0.11244 

JOINT 

302 

u- 

0.02742 

V« 

0.00824 

w> 

0.08451 

JOINT 

303 

u* 

-0.01478 

V« 

0.01766 

w« 

0.08521 

JOINT 

304 

u« 

-0.00324 

v» 

0.01437 

w> 

0.06831 

JOINT 

305 

u- 

-0.02757 

v« 

0. 

w» 

0.00149 

JOINT 

306 

u- 

0. 

V- 

0. 

w« 

0. 

REACTIONS  FOR  LOADING  CONDITION 

1  FOLLOW 

JOINT 

106 

RX- 

0. 

RY-  0.76800 

RZ- 

0.00000 

JOINT 

305 

RX» 

0. 

RY-  0.76800 

RZ- 

0. 

JOINT 

306 

RX« 

0.00000 

RY-  -0.00000 

RZ- 

0.00000 

REACTIONS  FOR  LOADING  CONDITION  2  FOLLOW 

JOINT 

106 

RX- 

0. 

RY-  -3.56212 

RZ> 

-2.91667 

JOINT 

305 

RX- 

0. 

RY-  0.76795 

RZ- 

0. 

JOINT 

306 

RX«  - 

10.00000 

RY-  4.33018 

RZ- 

2.91666 

REACTIONS  FOR  LOADING  CONDITION  3  FOLLOW 

JOINT 

106 

RX- 

0. 

RY-  3.16988 

RZ- 

4.33339 

JOINT 

305 

RX« 

0. 

RY-  46.37354 

RZ- 

0. 

JOINT 

306 

RX* 

-9.99999 

RY-  -34.54341 

RZ- 

2.16672 

REACTIONS  FOR  LOADING  CONDITION  4  FOLLOW 

JOINT 

106 

RX« 

0. 

RY-  -1.14434 

RZ- 

-1.41667 

JOINT 

305 

RX« 

0. 

RY-  -3.59811 

RZ- 

0. 

JOINT 

306 

RX« 

-1.00000 

RY-  2.74245 

RZ- 

-1.58334 
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TINSION  «VC  CONMESSION-¥E/ 

IM  LOAOINE  CONDITIONS- 


NUMtEN 

1 

2 

3 

4 

lit 

112 

-l.lMtlA 

-••fMll7 

-■.Iflfll 

-1.615111 

ifi 

113 

-1.116116 

-1.116116 

-1.111115 

1.314911 

111 

114 

-1.114132 

1.611179 

-1.396136 

1.582557 

ifi 

314 

-1.114132 

1.611179 

-1.396136 

1.249224 

lU 

113 

-1.114343 

-1.129684 

1.941263 

-1.182114 

if2 

114 

1.116111 

-1.631796 

-1.691678 

1.593117 

IfE 

115 

-1.141121 

1.131211 

-2.976561 

1.261853 

lf2 

116 

1.123157 

-1.558158 

-1.614582 

1.232224 

112 

314 

1.167164 

1.715667 

1.691678 

1.515778 

Ifl 

316 

1.134111 

-2.649821 

-6.253517 

1.544996 

IfS 

114 

1.167164 

-1.176774 

2.135658 

-1.347429 

IfS 

IIS 

1.134112 

1.388948 

3.541526 

-1.612211 

IfS 

314 

1.116111 

1.161647 

1.433783 

-1.149235 

IfS 

315 

1.123158 

1.161919 

9.751881 

-1.125523 

its 

316 

-1.141121 

1.131211 

-6.119461 

1.261852 

114 

115 

-1.115331 

-1.782259 

-1.292171 

1.655824 

lf4 

116 

1.113762 

1.543747 

2.217724 

-1.815669 

if4 

211 

1.142451 

1.224311 

1.743591 

-1.175862 

lt4 

212 

1.137141 

-1.417318 

-1.357147 

-1.383816 

211 

214 

1.152682 

1.234527 

1.743583 

-1.175861 

2i2 

213 

1.131149 

1.617853 

1.467465 

-1.118711 

2f2 

214 

-1.179585 

-1.639735 

-2.647482 

8.383816 

2i2 

215 

-1.166673 

5.456656 

1.156913 

1.159138 

2f3 

214 

-1.147622 

-1.417114 

1.584684 

-1.185246 

2t3 

215 

1.146463 

-3.684478 

-1.658697 

-1.165825 

2f4 

215 

1.111861 

-2.719428 

1.111614 

1.197816 

214 

216 

1.111161 

1.719926 

1.924447 

-1.256117 

2t4 

311 

1.852682 

-2.446775 

-1.942181 

-1.241978 

2i4 

312 

-1.147622 

-1.121688 

1.935216 

1.161529 

2t4 

313 

-1.179585 

1.186129 

-1.384123 

1.176761 

2tS 

216 

1.121145 

-1.971314 

-1.646411 

-1.173564 

2tS 

313 

-1.816788 

1.492156 

-1.561141 

-1.109723 

216 

312 

8.146464 

-2.731582 

6.194549 

-1.291697 

216 

313 

-1.866673 

2.441148 

-7.549448 

0.457171 

311 

312 

8.813418 

1.813419 

-1.111121 

0.110001 

311 

313 

1.883418 

1.113419 

1.111117 

-1.000000 

311 

314 

1.842458 

-2.457118 

-1.942178 

-0.241978 

312 

313 

8.831149 

1.985889 

1.221653 

-0.144568 

114 

213 

1.885186 

8.364668 

-1.584685 

0.085247 

115 

116 

8.817267 

8.287416 

1.299628 

-0.264698 

115 

213 

1.835578 

1.466539 

1.716713 

-0.363786 

116 

213 

-8.866672 

-5.884276 

-1.478284 

0.159135 

116 

212 

-1.217365 

5.378732 

-1.142281 

1.836816 

211 

212 

1.883418 

8.883416 

-1.118111 

0.000001 

211 

213 

1.883488 

1.883482 

1.111115 

0.000001 

312 

314 

8.815186 

8.978252 

-1.935214 

-0.161530 

312 

315 

-1.866673 

-2.831315 

-3.861149 

0.457172 

312 

316 

1.835581 

-1.738923 

11.312914 

-0.456607 

313 

314 

1.837149 

1.776161 

4.388749 

-0.176761 

313 

315 

-1.217366 

2.452536 

-13.359841 

1.618135 

314 

315 

8.813762 

-8.154852 

-1.491415 

0.601445 

314 

316 

-1.815331 

1.189837 

3.574659 

-1.179792 

315 

316 

1.817267 

8.911238 

1.634329 

-1.311114 

212  216 

13  tIARSCL 

14  fCHECR 

-1.816788 

-8.415652 

-1.167491 

0.142012 
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1MSI4  STAIR  MOtLfR  NUMBER  2lf 
SUMS  Of  THE  INTERNAL  STRESSES  AT  THE  JOINTS  FOLLOW 
LOADING  CONDITION  1 


JOINT 

FX 

FT 

FI 

111 

B.lltff 

-0.13542 

-0.01010 

112 

•  .lltlB 

-0.10461 

-0.00000 

113 

l.flBIB 

-0.10461 

0.00000 

114 

-B.ffifl 

-0.09787 

0.00000 

IfS 

't.flllf 

-0.08475 

-0.00000 

116 

f.ilflf 

0.68325 

0.00000 

211 

-C.lffll 

-0.03542 

0.00000 

212 

f.llfll 

-0.10460 

-0.00000 

213 

C.RIIfl 

-0.10460 

-0.00000 

214 

-t.fllfl 

-0.09787 

0.00000 

215 

-f.ffiff 

-0.08475 

-0.00000 

216 

-C.IIISB 

-0.08475 

-0.00000 

311 

B.lIBBB 

-0.03542 

0.00000 

3S2 

-B.llifl 

-0.10461 

-0.00000 

313 

e.iicie 

-0.10461 

0.00000 

314 

-0.09787 

-0.00000 

315 

-C.IICfl 

0.68325 

0.00000 

316 

C.eiBAI 

-0.08476 

0.00000 

LOADING  CCNOITICN  2 


JOINT 

FX 

FV 

FI 

101 

-0.00000 

-0.03542 

0.00000 

102 

-0.00000 

-8.10460 

0.00000 

103 

-0.00000 

-0.18461 

-0.00000 

104 

0.00000 

-0.09786 

-0.00001 

105 

0.00000 

-0.08475 

0.00000 

106 

0.00000 

-3.64687 

-2.91667 

211 

0.00000 

-0.03543 

-0.00002 

212 

4.99996 

-0.10461 

7.50002 

203 

0.00001 

-0.10464 

0.00002 

204 

0.00002 

-0.09786 

-0.00000 

205 

£.00001 

-0.08474 

8.00000 

206 

0.00000 

-0.08475 

0.00000 

301 

-0.00000 

-8.03542 

-0.00000 

302 

-0.00000 

-0.10460 

-0.00000 

303 

5.00000 

-0.10460 

-7.50000 

304 

-0.00000 

-0.09786 

0.00000 

305 

0.00000 

0.68320 

-0.00000 

306 

-10.00000 

4. 24542 

2.91666 

LOADING  CCNDITICN  3 

JOINT 

FX 

FY 

FI 

101 

-0.00000 

-0.00002 

-0.00001 

102 

-0.00000 

0.00000 

-0.00003 

103 

-0.00000 

-4.99998 

0.00003 

104 

0.00000 

-0.00001 

0.00000 

115 

-0.00000 

-0.00000 

-0.00000 

116 

-0.00001 

3.16986 

4.33340 

201 

-0.00000 

0.00003 

0.00002 

212 

5.00004 

-0.00001 

-0.00008 

203 

0.00081 

0.00004 

0.00002 

204 

0.00000 

-10.00003 

0.00000 

205 

-0.80002 

-0.00002 

-0.00006 

206 

-0.00000 

0.00001 

-7.50000 

301 

0.00003 

0.08002 

0.00004 
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3i2 

-(.00000 

-0.00000 

-0.00014 

333 

5.00003 

-0.00002 

1.00006 

3f4 

0.00001 

-0.00001 

-0.00000 

3i5 

-0.00001 

46.37355 

0.00000 

3f6 

-5.49999 

-34.54342 

2.16672 

LOADING  CCNOITIOK  4 

JOINT 

FX 

FY 

FZ 

101 

1.00000 

2.00000 

3.00000 

102 

0.00000 

-0.00000 

0.00800 

103 

0.00000 

-0.00000 

-0.00000 

104 

0.00000 

0.00000 

0.00000 

105 

0.00000 

0.00000 

0.00000 

106 

(.00000 

-1.14434 

-1.41667 

201 

0.00000 

-0.00000 

-0.00000 

202 

-0.00001 

-0.00000 

0.00001 

203 

-0.00000 

-0.00000 

-0.00000 

204 

(.00000 

0.(0000 

0.00000 

205 

0.00000 

0.00000 

0.00000 

206 

0.00000 

-0.00000 

-0.00000 

301 

-0.00000 

-0.00000 

-0.00000 

302 

(.00000 

0.00000 

0.00001 

303 

-0.00000 

0.00000 

-0.00001 

304 

-0.00000 

0.00000 

0.00000 

305 

0.00000 

-3.59811 

-0.00000 

306 

-1.00000 

2.74245 

-1.58334 
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CHAPTER  3 
OPERATOR'S  MANUAL 

A.  INTRODUCTION 

Aspects  of  the  STAIR  System  which  affect  the  computer  operator  and  programmer  are  dis¬ 
cussed  in  this  and  the  following  chapter.  The  principal  item  covered  here  is  running  a  STAIR 
problem  from  data  sheets  prepared  by  the  user. 

The  machine  procedure  is  described  for  the  IBM-709  and  -7090.  A  similar  procedure  is 
used  on  the  IBM-704  except  for  tape  assignments. 

B.  INPUT-OUTPUT 

Data  prepared  by  the  user  is  punched  into  two  card  decks  (structural  data  for  STAIR  and  load 
data  for  SLOP).  Both  decks  are  converted  directly  to  BCD  (Binary  Coded  Decimal)  input  tapes. 
The  load  data  program  binary  must  first  be  run  on  line  with  a  load  assembly  program  (SLOP) 
tape  (B5),  which  assembles  the  loads  on  a  binary  input  tape.  All  input  to  the  STAIR  program  is 
thus  from  tape. 

Off-line  output  from  tape  B2  is  printed  on  peripheral  equipment  under  program  control* 
after  a  successful  run.  On-line  output  includes  error  diagnostics  and  tracking  information,  de¬ 
pending  on  the  sense  switch  settings. 

C.  MACHINE  OPERATION  | 

Starting  from  prepared  data  and  program  tapes,  a  STAIR  problem  is  run  as  follows; 

1.  Mount  tapes. 

B2  Output  (any  tape) 

A6  Buffer  (any  tape) 

A  7  Buffer  (any  tape) 

A8  Load  Data  (binary-output  from  SLOP) 

A9  Structural  Data  (cards  pre-stored  off-line) 

B4  STAIR  Program  Tape 

2.  Run  program  by  loading  STAIR  loader  on  line  (one-card  tape  loader, 

T/4  transfer  card  and  several  blank  cards). 

3.  Dismount  and  save  structural  data,  load  data  and  program  tapes. 

4.  List  off-line  output. 

D.  SENSE  SWITCHES  AND  INDICATOR  LIGHTS 

Sense  switch  settings  allow  the  following  options  during  the  running  of  a  problem; 

SS  1  UP  Stop  at  all  detected  errors  after  printing  diagnostic. 

DN  Print  diagnostic  but  do  not  stop  after  certain  errors  are 

detected  in  SLOP  and  subroutine  STRIX  (see  User's  Manual). 


"■Experience  has  shown  that  it  is  more  convenient  to  have  the  output  printed  single-spaced,  but 
this  varies  with  the  peripheral  equipment  available. 


SS  2  UP 
DN 

SS  3  UP 
DN 

SS  6  UP 
DN 


No  tracing  print  out. 

Print  tracing  information  after  completion  of  each  step  in 
ANALYSIS  PROGRAM. 

No  off-line  output  from  SLOP. 

Write  assembled  loads  on  output  tape  B2  for  listing. 

No  off-line  output  from  intermediate  steps. 

Write  all  information  completed  in  this  step  on  output  tape 
B2  for  listing  (used  only  for  program  degugging). 


Under  no  circumstances  should  the  I-O  check  light  be  on  at  the  completion  of  a  successful 
run.  The  Divide-Check  and  Overflow  indicators  may  be  on. 


E.  PROGRAM  STOPS 

Most  program  stops  will  occur  after  an  on-line  diagnostic  is  printed.  If  a  stop  occurs  with¬ 
out  a  diagnostic,  the  Instruction  Location  Counter  (ILC)  should  be  recorded  and  a  programmer 
consulted.* 

At  the  end  of  a  successful  run,  all  tapes  on  DSC  A  are  rewound,  and  the  program  is  in  po¬ 
sition  to  start  a  new  problem.'  The  next  problem  may  be  run  by  mounting  new  data  tapes  AS 
and  A9,  rewinding  tape  B4,  replacing  the  STAIR  loader  in  the  card  reader  and  then  clearing  and 
loading  cards. 

After  stops  caused  by  tape  failure,  the  program  may  be  rerun,  but  a  different  tape  and  tape 
drive  should  be  tried.  If  repeated  stops  occur,  one  of  the  error  stops  without  a  diagnostic  might 
be  suspected. 

Errors  in  data  input  must  be  corrected  in  the  appropriate  card  deck,  and  that  deck  must  be 
rewritten  on  tape  before  the  run  is  repeated. 


P.  INTERRUPTION  OF  A  RUN 

For  problems  requiring  a  large  amount  of  machine  time  for  their  solution,  it  may  not  be 
possible  to  do  all  the  computations  without  an  interruption  which  would  require  getting  off  the 
machine. 

To  interrupt  the  program: 

1.  Mount  a  tape  on  B5  for  dumping  the  core. 

2.  Put  Sense  Switch  4  down.  This  will  cause  the  master  routine  to  take 

a  core  dump  on  tape  BS  after  the  current  step  is  completed,  then  stop. 

3.  The  partial  output  is  on  tape  B2.  Save  tapes  B5,  A6,  A7,  A8,  A9  and  B4. 

To  continue  the  run  at  a  later  time; 

1.  Rewrite  the  structural  data  tape  A9  (off-line)  with  the  data  for  the 
remaining  steps. 

2.  Mount  tapes  BS,  A6,  A7,  A8,  A9,  B4  and  a  new  B2. 

3.  Put  Sense  Switch  5  down. 

4.  Run  program  by  loading  STAIR  loader  (tape  B4  card  followed  by  a 
T/4  card  and  several  blank  cards). 


*Also  see  p.  149. 


CHAPTER  4 

PROGRAMMER'S  MANUAL 


A.  INTRODUCTION 

The  STAIR  System  consists  of  11  routines,  written  on  a  program  tape  and  obtained  by  the 
computer  during  an  analysis  problem,  and  a  load  assembly  program  used  to  prepare  data  for 
the  analysis. 

The  routines  on  the  program  tape  may  be  grouped  as  follows; 

1.  STAIR  operating  routines  —  nine  routines  (plus  a  subroutine)  which  per¬ 
form  the  operations  described  in  Chapter  1. 

2.  STAIR  service  routines  -  two  routines  which  are  called  by  several  of  the 
operating  routines,  and  a  master  routine  which  initiates  each  analysis 
problem  and  calls  successive  operating  routines. 

3.  Library  routines  —  conversion,  I-O  and  elementary  function  routines,  as 
punched  from  the  FORTRAN  II  compiler. 

4.  MISLAM  —  an  executive  routine  which  writes  the  STAIR  program  tape 
from  a  binary  deck  and  controls  its  operation. 

Except  for  the  library  routines,  all  these  routines  and  a  pair  of  load  assembly  programs 
are  described  on  the  following  pages. 

Since  the  number  of  instructions  comprising  the  STAIR  routines  is  about  30,000,  there  is 
not  enough  storage  space  in  a  32,767-register  computer  tc  store  all  the  programs  and  a  large 
matrix,  simultaneously.  Instead,  the  nine  operating  routines  are  stored  on  tape  and  called  into 
core  only  when  needed. 

Transfer  of  control  between  operations  is  handled  by  the  master  routine.  Common  storage 
(upper  memory)  is  used  extensively  to  allow  continuity  between  routines. 

B.  DEFINITION  OF  MATRIX  TERMS 

A  stiffness  matrix  K  is  formed  from  data  describing  the  geometry  of  a  structure  and  modu¬ 
lus  of  elasticity  of  the  material  E. 

A  matrix  equation  is  composed  of  a  stiffness  matrix  and  a  load  matrix  F  3;thich  determine 
a  displacement  matrix  U  =  K*^F  or  F  =  KU.  K  is  a  symmetrical  matrix. 

Where  displacements  are  governed  by  external  restraints,  rows  corresponding  to  those 
restrained  displacements  may  be  removed  from  the  stiffness  and  load  matrices  and  stored  as 
reaction  matrices  (Kp,  Fp).  Reactions  may  be  evaluated  from: 

R  =  KpU4Fp  . 

A  matrix  equation  may  be  partitioned  into  the  following  form: 


1 

iT 

_ 

- 1 

c 

_ 1 

- 1 

_ 1 

•^21  i  K22 

^2 

/2 

Block  elimination  yields  an  explicit  expression  for  U2: 


^!z-^zz-^Iz^u  ^iZ 


f|=F3-kT 


-T  . 


Kj2  transpose  of 


The  superscript  e  denotes  an  "effective"  matrix.  The  eliminated  displacements  may  be 
evaluated  from  a  group  of  matrices  called  "residuals"; 


U,  =  K-  (F,  -  K,2  U^) 


C.  CORE  STORAGE  ALLOCATION 

A  map  of  core  storage  during  an  analysis  problem  is  given,  with  approximate  actual  locations. 


Buffer  Storage 


Common 

Storage 

(Length  Voriet) 


Not  Used 


STAIR  Operating  Routine 
(Length  Varies) 


STAIR  Service  Routines 


Library  Routines 


MASTER  Routine 


MISLAM 


Loader  BSS 


Location  (octal) 


77777 

77461 


11 

0 


Common  storage  contains  a  list  of  parameters  and  arrays  to  be  defined  in  the  next  section. 
All  arrays  are  FORTRAN-stored,  one-dimensional  arrays.  Only  the  upper  half  (above  the  main 
diagonal)  of  a  symmetric  matrix  is  stored,  so  addresses  are  computed  for  a  triangular  array. 
Stiffness  matrices  are  stored  row-wise  and  load  matrices  are  stored  column-wise. 
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D.  DEFINITION  OF  COMMON  VARIABLES 


Common  storage  contains  the  following  list  of  variables  during  execution  of  the  operating 
routines; 


IDP 

NSTEP 

E 

NREF 

NJREAC 

NUT 

NK 

NU 

NJT 

NJSEC 

NRJSEC 

NRCSEC 

JREAC 


NTAPE 

NROW 

STRIX 


problem  identification  number 

consecutive  number  of  steps  in  ANALYSIS  PROGRAM  being 
executed 

elastic  modulus  of  bar  material 

number  of  residual  and  reaction  matrices  written  on  tape  A6 
number  of  reaction  joints  in  structure 
number  of  units  in  structure 
number  of  loading  conditions 

number  of  unit  whose  matrix  is  in  primary  storage  area 
number  of  joints  in  primary  row  sequence 

number  of  joints  being  eliminated  from  primary  row  sequence, 
or  number  of  joints  in  secondary  row  sequence 

number  of  reaction  joints  among  joints  being  eliminated 
number  of  reaction  components  among  joints  being  eliminated 

list  of  all  reactions  acting  on  structure  (four  words  per  reaction 
joint  giving  joint  number  and  code  indicating  whether  or  not  com¬ 
ponents  of  reaction  act  in  x-,  y-  or  z-directions);  (0)  means  no 
reaction,  (1)  means  a  reaction 

tape  reading  error  indicator  (set  by  FIXTPE  and  sensed  by 
routine  calling  FIXTPE) 

list  of  joint  numbers  associated  with  matrix  in  primary  storage 
area 

primary  matrix  storage  area 


The  following  arrays  are  included  in  common  storage  by  all  routines  except  STRIX  (to  be  dis¬ 
tinguished  from  array  STRIX),  BARSOL  and  CHECK. 


NROSEC  secondary  row  sequence,  associated  with  matrix  in  secondary 
storage  area 

LRSEC  list  of  reactions  among  joints  being  eliminated 
STSEC  secondary  matrix  storage  area 


Additional  parameters  and  arrays  may  be  included  in  common  storage  by  the  use  of  certain  op¬ 
erating  routines,  but  these  parameters  do  not  affect  the  continuity  of  operation  among  routines. 


E.  TAPE  STORAGE  FORMATS 

Storage  formats  for  the  input  and  buffer  tapes  attached  to  DSC  A  are  given  below.  The 
format  of  the  program  tape  B4  will  be  described  with  MISLAM. 

Tape  A6  —  Residuals  and  reaction  matrices  in  two  binary-record  groups. 

Record  1  -  IDOUT,  NJT,  NJSEC,  NRJSEC,  NRCSEC 
Record  2  -  NROW,  LRSEC,  STSEC 

Tape  A7  —  Matrix  equation  storage  in  three  binary-record  groups. 

Record  1  -  NHED,  NJT.  NJSEC,  NRJSEC,  NRCSEC,  NU 

Record  2  -  NROW,  STRDC 

Record  3  -  NROSEC,  LRSEC,  STSEC 
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Tape  A8  —  Load  data  with  consecutive  binary  records  containing  the  load 

matrix  for  consecutive  units  of  the  structure.  The  first  word  of 
each  record  is  a  heading  NHBD. 

Tape  A9  —  Structural  data  in  groups  of  BCD  records. 

Following  an  initial  group  of  two  data  records,  the  records  for 
the  step  number  and  the  c^eration  number  follow,  each  of  which 
is  followed  by  the  data  records  for  the  step. 

Records  1  and  2  —  General  structure  data 
Record  3  —  Step  number,  operation  munber 
Record  4. ...  n  —  Data  for  this  step 
Record  n+1  -  Step  number,  operation  number 

NHED  identifies  a  record  group  u  as  32767-u. 

IDOUT  identifies  a  two-record  group  as  a  residual  (0)  or  reaction  matrix  (1). 

F.  STAIR  OPERATING  ROUTINES 

The  nine  operating  routines  are  described  on  the  following  pages.  Each  description  has  the 
following  format: 

(a)  Operation  —  definition  of  the  matrix  operation  performed. 

(b)  Mathematical  method  —  brief  description  of  technique  for  performing 
the  operation. 

(c)  Initial  condition  of  core  —  common  variables  which  must  be  present  to 
maintain  continuity. 

(d)  Error  stops. 

(e)  Output  -  to  tape  only;  on-line  output  consists  of  error  diagnostics  and 
tracing  information. 

(f )  Final  condition  of  core  -  common  variables  which  are  changed,  affecting 
continuity. 

(g)  Addressing  diagrams. 

(h)  Block  diagrams  —  flow  charts. 


1.  STRDC 

Operation:  Forms  a  stiffness  matrix  equation  in  core  from  input  defining  the  geometry  of 
a  structure  and  its  load  matrix. 

Mathematical  Method:  Terms  of  the  matrices  for  single  bars  are  computed  and  superim¬ 
posed  on  the  array  of  the  final  matrix.  The  location  of  the  final  matrix  array  of  the  bar  con¬ 
necting  joints  a  and  b  is  determined  by  the  positions  of  a  and  b  in  the  list  of  joint  numbers 
NROW.  The  load  matrix  for  the  structure  is  read  directly  from  tape  A8,  and  is  stored  imme¬ 
diately  after  the  stiffness  matrix.  A  vector  check*  indicates  if  all  the  bars  meeting  at  any  joint 
lie  in  a  plane;  such  a  geometry  would  result  in  an  ill-conditioned  stiffness  matrix. 

Initial  Condition  of  Core:  Common  variables  E,  NUT  and  NK  must  be  predetermined. 


*The  check  for  a  planar  joint  is  performed  in  the  following  manner:  Let  =  unit  vector  in  di¬ 
rection  of  bar  1;  L2  =  unit  vector  in  direction  of  bar  2;  L  =  unit  vector  in  direction  of  any  other 
bar.  First,  the  vector  product  C  -  Lj  x  L2  is  executed.  If  each  component  of  C  is  less  than 
0.000001,  the  bars  are  considered  collinear,  and  the  unit  vector  associated  with  the  next  bar  is 
used  as  112  to  form  a  newC.  Then,  the  scalar  product  /C-L/  is  executed.  If  /C-L/ <  0.00001 
for  every  remaining  bar,  the  joint  is  considered  planar. 


r 


I 


Error  Stops: 

(a)  Number  of  terms  in  matrix  equation  >  16000. 

(b)  Number  of  bars  in  unit  structure  >  600. 

(c)  Unit  number  >  number  of  units. 

(d)  A  bar  number  includes  the  number  of  a  joint  now  in  NROW. 

(e)  All  the  bars  meeting  at  a  joint  lie  in  a  plane  (see  Sense  Switch  1  settings). 

(f )  Tape  failure  searching  for  load  matrix. 

Output:  None. 

Final  Condition  of  Core:  Common  variables  STRK,  NROW,  NJT  and  NU  define  the  mat."ix 
equation  which  has  been  formed.  Common  variables  NJSEC,  NRJSEC  and  NRCSEC  are  zeroed. 


i 


(? 
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2.  EFFRIX 


Operation;  Eliminates  a  block  of  rows  from  a  matrix  equation  in  core. 

Mathematical  Method:  The  rows  to  be  eliminated  must  be  at  the  top  of  the  original  matrix 
equation,  so  that  it  may  be  subdivided  as  indicated: 


'Kll 

_ 1 

P 

Ha 

^21 

^22 

^2 

^2 

. 

. 

If  any  of  the  displacements  U  are  fixed  by  external  reactions,  the  corresponding  rows  from 

and  are  then  separated,  condensed  and  re-stored  with  the  reaction  rows  removed, 
is  expanded  into  square  form  in  the  area  previously  occupied  by  and  inverted. 

Effective  matrices  are  formed  according  to  the  following  equations  and  stored  at  the  head  of  the 
array  occupied  by  the  original  matrix  equation. 


k|2  =  K22-kT  K^-^K 


-i 


Finally,  the  residuals 


12 


are  written  on  tape  A7. 


Initial  Condition  of  Core:  Common  variables  JREAC,  NJREAC  and  NK  must  be  present.  A 
matrix  equation,  as  formed  by  STRDC  or  MADD,  must  be  in  array  STRIX,  with  identifying 
parameters  in  NROW  and  NJT. 


Error  Stops:  Violations  of  size  restrictions,  as  indicated  in  the  addressing  diagrams. 

Output:  A  reaction  matrix,  if  any  of  the  eliminated  displacements  correspond  to  reactions 
and  residuals  (all  on  tape  A6). 


Final  Condition  of  Core:  Effective  matrices  are  stored  in  STRDC,  and  residuals  in  STSEC. 
The  effective  matrix  equation  is  identified  by  NJT  and  NROW,  both  of  which  are  modified  by  the 
elimination  of  joints.  NROSEC  contains  the  original  NROW;  NJSEC  contains  the  number  of  joints 
eliminated  from  NROW.  The  location  and  number  of  reactions  among  the  eliminated  displace¬ 
ments  are  indicated  by  NRJSEC,  NRCSEC  and  LRSEC. 
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NRA«3NJSEC-NRCSEC 

NRB=3(NJT-NJSEC) 

SIZE  RESTRICTIONS; 

NRA  (NRA  +  l)/2  +  NRA  (NRB  +  NK)  ^  9500 

NRA*  ^  3  NJSEC  (3  NJSEC  +  l)/2  +  9  NJSEC  (NJT  -  NJSEC) 
(3  NJT  +  NRB  +  l)/2  3  NJSEC  >  (NRB  +  IV'2  NRB 

MATRIX  RESIDUAL  STORAGE  BY  EFFRIX 
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1.  CHECK  JOINTS  TO  BE  ELIMINATED 
FOR  REACTION  COMPONENTS 


CONTINUE 


18 


».  FORM  EFFECTIVE 
LOAD  MATRIX  F^ 


LSLS 

-  NRPX(NRP+  11/2 

+  NRS 

LBLS 

-  LBS  +  NRA  X  NRB 

LS- 

LSLS  +  1 

LT- 

LT+  1 

SET  K- 


STEP  K  +  1 


]-| 


SET  JDEF  -  I 


SET  JDEF  =  1 


SET  XS  =■  1  I 

”1_ 


IS  -  IS  +  NRS 

I 


TERMINATE  ON 
K-NK 

YES  NO 
□ - 


STEP  JDEF  +  L 


TERMINATE  ON 
JDEF  •  NRA 
YES  NO 


STEP  JDEF  +  1 


LS  °  LBLS  + 

(K-  l)XNRA+  1 
LA  °  JDEF 

zn 


SET  JRB  =  1 


TERMINATE  ON 
JDEF  >  NRA 
YES  NO 

1 - c 


STEP  XS  +  1 


Lt  =  LBS  +  XS 
STRIX  (LT)  - 
STRIX  (LS) 

~T — : 


SET  JRB  -  1 


LT  -  LT  +  I 
LS  »  LS  +  1 

I 


TERMINATE  ON 

XS  -  NRB 

YES  NO 


DEFSEC  (JDEF)  -  0 

STEP  JRB  +  1 

DEFSEC(XEF)  -  DEFSEC(JDEF) 
+  STSEC(LA)  XSTSEC(LB) 

~1 - 

I  JRB- JDEF 


LA  -  LA  +  NRA 
-  JRB 

_ 1 

LA-LAH 

- 1 

I  LB  =  L  B  +  1~| 


TERMINATE  ON 
JRB  >  NRA 
YES  NO 


STEP  JRB  + 1 


STRIX  (LT)  -  STRIX  (LT) 

-  DEFSEC(JRB)  X  STSEC(LB) 


I 


LB  =  LB  +  NRB 


TERMINATE  ON 
JRB -NRA 
YES  NO 


I  CONTINUE  I 
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3.  MATRO 


f 


I 


f 


Operation:  Forms  the  stiffness  matrix  equation  of  a  unit  from  the  stiffness  matrix  of 
another  geometrically  similar  unit  and  the  load  matrix  of  the  new  unit. 

Mathematical  Method:  Transformation  formulae  are  applied  to  3  X  3  blocks  of  terms  to  ac¬ 
complish  a  rotation  in  the  x-y  coordinate  plane.  The  load  matrix  of  the  new  unit  is  read  directly 
from  tape  A8.  If  an  effective  matrix  is  rotated,  its  residuals  are  also  rotated  and  a  new  effec¬ 
tive  load  matrix  is  computed.  Since  the  transformation  formulae  are  applied  to  3  X  3  blocks, 
they  may  not  be  applied  to  a  residual  from  which  rows  have  been  removed  to  form  a  reaction 
matrix. 


Initial  Condition  of  Core:  The  matrix  equation  of  a  single  unit  of  a  structure  must  be  stored 
in  core,  as  by  STRDC  or  EFFRDC. 

Error  Stops: 

(a)  New  unit  number  >  quantity  of  units. 

(b)  Reaction  joints  among  joints  in  old  or  new  effective  matrices. 

(c)  Tape  failure  searching  for  load  matrix. 

Output:  Residuals  for  the  new  unit  are  written  on  tape  A6  when  an  effective  matrix  is 
transformed. 


Final  Condition  of  Core:  The  general  map  of  core  storage  is  unchanged;  transformed 
matrices  occupy  STRDC  and  STSEC,  and  the  new  unit  is  identified  by  NU,  NROW  and  NROSEC. 

Addressing  Diagram:  Terms  of  a  3  x  3  block  are  addressed  from  the  following  location 
symbols  computed  on  the  diagonal  of  the  block: 


LA 

LD 

LF 

% 


AAATRO-GENERAL 
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TERMINATE  ON 
K-NK 

YES  NO 


TERMINATE  ON 
NR  =  NRP 
YES  NO 
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4.  MADD 


Operation:  Superimposes  the  stiffness  matrix  equations  on  two  adjacent  parts  of  a  structure. 

Mathematical  Method:  A  matrix  equation  initially  in  STRK  is  transferred  to  STSEC,  and 
STRDC  is  cleared.  The  matrix  in  STSEC  is  identified  by  NJSEC  and  NROSEC.  A  new  list  of  NJT 
joint  numbers  is  read  into  NROW  which  identifies  the  new  matrix  equation  to  be  formed  in  STRDC. 
Joint  numbers  in  NROSEC  are  compared  to  those  in  NROW;  the  consecutive  number  of  a  joint  in 
NROW  replaces  the  corresponding  joint  number  in  NROSEC  and  3X3  blocks  of  terms  in  STSEC 
are  superimposed  on  the  STRDC  array,  with  addresses  computed  from  the  renumbered  NROSEC. 
NK  X  3  blocks  from  the  load  matrix  in  STSEC  are  stored  similarly.  When  the  original  in-core 
matrix  has  been  restored,  a  second  matrix  is  located  on  tape,  read  into  STSEC,  NROSEC  and 
NJSEC,  and  the  superposition  process  is  repeated. 

Initial  Condition  of  Core:  A  matrix  equation  as  formed  by  STRDC,  MATRO,  EFFRDC  or 
another  MADD  must  be  in  storage  area  STRDC;  it  is  identified  by  NROW  and  NJT. 

Error  Stops: 

(a)  Number  of  registers  occupied  by  final  matrix  equation  >  16000. 

(b)  Number  of  registers  occupied  by  either  equation  being  added  >  9500. 

(c)  A  joint  number  in  NROSEC  which  is  not  in  NROW. 

(d)  Tape  failure  searching  for  matrix  on  tape,  or  no  matrix  on  tape 
corresponding  to  the  tape  location  given  as  input. 

Output:  None. 

Final  Condition  of  Core:  A  matrix  equation  is  stored  in  STRDC,  with  identifying  parameters 
in  NROW  and  NJT.  (Common  assignments  for  STSEC  and  NROSEC  differ  in  their  routine  from 
the  others,  but  they  do  not  affect  continuity.)  NJSEC  is  zeroed. 


MAOO  ADDRESSING  DIAGRAM 


MATRIX  SUPERPOSITION  BY  MADD 


READ  NJTF,  MATLOC 


MAT-0 


INITIALIZE  STORAGE  OF 
IN-CORE  MATRIX 


READ  TAPE  7 
NHED,  NJSEC 


2.  RENUMBER  SECONDARY  ROW 
SEQUENCE  (NROSEC) 


NRS  -  3NJSEC 
MAXS-NRSX(NRS  +  l)/2 
+  NRSXNK 


3.  SUPERIMPOSE  SECONDARY 
MATRIX  ON  PRIMARY 
MATRIX 


9500  -  MAXS 


4.  SUPERIMPOSE  SECONDARY 
LOAD  MATRIX  ON  PRIMARY 
LOAD  MATRIX 


READ  TAPE  7 
NROSEC  (I),  STSEC  (1) 


r±^ 

1  MAT  1 

1  i 

1  W 

NJSEC -0 
NRJSEC  -  0 
NRCSEC-0 


MAT  -  MAT +1 


PRINT  BtROR: 

SIZE  LIMITATIONS 
EXCEEDED 


3.  SUPERIMPOSE  SECONDARY 
MATRIX  ON  PRIMARY 
MATRIX 


STEP  I  +  I 


SET  J  - 1 


STEP  J  +  I 


NROSEC(l)-NROSEC(J) 


>-NROSEC  (I) 
'  -  NROSEC  (J) 


MP  -  NROSEC  (J) 
NP  -  NROSEC  (I) 


MP=NP-NR0SEC(I) 


K  =  (2NRS-3I  +  4)(I-I)V2 
+  3(J-I)+  1 

KT=K  +  3(NJSEC-l  +  1) 

KTT  =  KT+3(NJSEC-I+  1)-1 
L  =  (2NRP  -  3MP  +  4)  (MP  -  1  )3/2 
+  3(NP-MP)+  1 
LT=L  +  3(NJT-MP  +  1) 

Ln  =  LT+3{NJT-MP+1)-l 


NROSEC(I)  -  NROSEC(J) 


STRIX(L)  =  STRIX(l)+  STSEC(K) 

STRIX(LT)  =  STRIX(LT)  +  STSEC(KT) 

STRIX(LTT)=  STRIX(Ln)+  STSEC(KTT) 

STRIX(L  +  1 )  =  STRIXIL  +  1 )  +  STSEC(K  T  -  1 ) 
STRIXCL  +  2)  =  STRIX(L  +  2)  +  STSEC(Kn  -  2) 
STRIX  LT  - 1 )  =  STRIXILT  -  1 )  +  STSEC(K  +  1 ) 
STRIXILT  +  t)  =  STRIX(LT  +  1)  +  STSEC(Kn  -  1 ) 
STRIxan  -  2)  =  STRIXILTT  -  2)  +  STSECIK  +  2) 
STRIX(Ln-l)  =  STRIX(LTT-I)+  STSEC(KT+  1) 


STRIX(LT  -  1 ) «  STRIXdT  -  1 )  +  STSEC(KT  -  1 ) 
STRIXI  in  -  1 )  •  STRIXlin  -  1 )  +  STSECIKH  -  1 
STRIxan  -  2)  ■  STRIX(LTT  -  2)  +  STSEC(Kn  -  2 


STSEC(Kn- 


L)  -  STRIX(L)  +  STSEC(K) 

L+  1)-STRIX(L+  1)+  STSEC(K+  n 
L  +  2) »  STRIX(L  +  2)  +  STSEC(K  +  2) 
LT)-STRIX(LT)+STSEC(KT) 

LT  +  1 )  =  STRIXdT  +  1 )  +  STSEC(KT  +  1 ) 
LTT)  -  STRIXdn)  +  STSEC(Kn) 


CONTINUE 


TERMINATE  ON 
l-NJSEC 


TERMINATE  ON 
J  -  NJSEC 
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5.  MATOUT 


Operation:  Writes  a  stiffness  matrix  on  tape  from  core. 

Mathematical  Method:  Matrices  and  identifying  parameters  are  written  as  three  record 
groups,  in  the  format  of  Sec.  E. 

Initial  Condition  of  Core:  Matrix  equations  may  exist  in  core  as  stored  by  STRDC,  MATRO, 
EFFRK.  MADD,  MATIN  or  another  MATOUT. 

Error  Stop:  Tape  failure  searching  for  tape  location  preceding  location  in  which  MATOUT 
is  to  write,  or  attempt  to  find  a  tape  location  at  which  nothing  has  been  written. 

Output:  None. 

Final  Condition  of  Core:  No  change  from  start  of  operation,  except  that  NTAPE  is  nonzero 


if  tape  search  has  failed. 


MATOUT-MATIN  ARRAY  DIMENSIONS: 


( 


STRIX 


NROW 


NO.  OF  TERMS 


MAXM  -  3  NJT  (3  NJT  +  \V2 
+  3  NJTXNK 


NJT 


NRS  =  3  NJSEC-  NRCSEC 
MAXM  =  NRS  (NRS  +  1)/2 

+  NRS  (3  NJT  +  NX) 


NRJSEC 


NROSEC 


^  NJSEC 

NJT  ^ 

1 _ 

1 

MAXR  =  4  NRJSEC  +  1 


MAXJ  »  NJSEC  +  NJT 
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6.  MATIN 


Operation:  Reads  a  matrix  into  core  from  tape. 

Mathematical  Method:  Matrices  and  identifying,  parameters  are  read  into  core  from  three 
record  groups  in  the  format  of  Sec.  E. 

Initial  Condition  of  Core:  All  matrices  and  identifying  parameters  are  reset  by  MATIN. 

Error  Stop:  Tape  failure  searching  for  tape  location,  or  attempt  to  find  record  which  has 
not  been  written. 

Output:  None. 

Final  Condition  of  Core:  Same  as  condition  when  the  record  group  was  written  by  MATOUT 
(in  regard  to  continuity  of  operation). 

Array  Dimensions:  See  MATOUT  description. 

-  f  . 


7.  DEFSOL 


Operation:  Evaluates  joint  displacements  and  reactions  from  a  final  stiffness  matrix  equation 
in  core  and  a  series  of  residuals  and  reaction  matrices  on  tape. 

Mathematical  Method:  A  set  of  displacements  is  determined  from  the  matrix  equation  in 
STRK: 

KU  =  F  . 

First,  if  any  of  the  displacement  u  is  restrained  by  reactions,  the  row  corresponding  to  the 
reaction  component  is  removed  from  K  and  F  and  written  on  tape  as  reaction  matrices 
(Kj^,  Fj^).  K  and  F  are  condensed  and  re-stored  in  STSEC;  K  is  expanded  into  square  form  in 
STRDC  and  inverted.  The  displacement  matrix  U  is  formed  by  multiplication  of  K’^  and  F  and 
is  stored  in  STSEC.  An  addressing  system  is  established,  and  STRIX  is  initialized  to  receive 
displacements  and  reactions.  After  the  displacement  matrix  U  is  re-stored  into  STRDC,  suc¬ 
cessive  groups  of  reaction  matrices  or  residuals  are  read  back.  Reactions  are  computed  by 
the  equation 

R  =  Kjj  U  and  Fj^  . 

Displacements,  previously  eliminated,  are  evaluated  with  other  displacements,  which  were 
eliminated  at  a  later  step: 

Ui  =  k;*(F,  -  U,)  . 

Initial  Condition  of  Core:  A  matrix  equation  must  be  in  STRDC  as  left  by  any  of  the  pre¬ 
ceding  routines. 

Error  Stops: 

(a)  Violations  of  size  restrictions,  as  indicated  in  the  storage  diagrams. 

(b)  Tape  failure  or  machine  failure  in  comparing  reaction  joint  numbers, 
read  back  with  a  reaction  matrix,  to  reaction  joint  numbers  stored  in 
core. 

Oitput:  Reaction  matrix  if  reactions  are  located  among  final  displacements.  Displace¬ 
ments  and  reactions  are  written  on  the  output  tape  B2  in  the  same  format  in  which  they  are 
stored  in  core  (one  record  per  joint,  giving  joint  number  and  u,  v  and  w  displacements). 

Final  Condition  of  Core:  Displacements  and  reactions  are  left  in  STRDC  as  they  were  stored. 
The  only  operation  which  may  logically  follow  DEFSOL  is  BARSOL. 
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OEFSOL  MATRIX  ADDRESSING  DIAGRAM 


REAaiON  EQUATIONS  (STSEC) 
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(a)  Displacement  storage  for  one  loading 

M.  =  number  of  joint  numbers  prefixed  by  unit  number  i. 


Displacements  for  loading  1 


Displacements  for  loading  NK 
Reactions  for  loading  1 


Reactions  for  loading  NK 


(b)  STRDC  map  for  NK  loadings 

4  X  NK  X  (number  of  joints  -f  number  of  reactions)  <16000. 


CONTINUE 


CONTINUE 


PRINT  ERROR:  SIZE 
LIMITATIONS  ARE 
EXCEEDED 


STOP 


CONTINUE 
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8.  BARSOL 


Operation:  Computes  the  stress  in  any  bar  from  the  displacements  and  coordinates  of  the 
joints  at  its  ends. 

Mathematical  Method:  Unit  stresses  are  computed  and  stored  for  all  loading  conditions  for 
consecutive  bars  in  an  input  list  of  bar  numbers. 

Initial  Condition  of  Core:  Common  variables  NK  and  NUT  must  be  present.  Joint  displace¬ 
ments  in  STRDC  and  the  addressing  array  lUADR  must  be  as  left  by  DEFSOL. 

Error  Stops: 

(a)  Number  of  bars  >1000. 

(b)  Number  of  bars  x  NK  >4300. 

(c)  Number  of  joints  at  ends  of  bars  >800. 

(d)  No  stop  occurs  if  an  illegal  bar  number  is  given  or  if  a  needed  set 

of  coordinates  is  not  supplied,  but  the  error  is  indicated  in  the  output. 

Output:  Bar  stresses  are  written  on  tape  B2  for  listing. 

Final  Condition  of  Core:  Common  locations  above  STRIX  are  not  changed. 
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CONTINUE 
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9.  CHECK 


Operation!  Sums  the  bar  forces  in  the  three  coordinate  directions  at  each  joint  for  each 
loading  condition. 

Mathematical  Method:  The  x,  y  and  z  forces  in  each  bar  are  read  from  tape  A7  and  stqier- 
impoaed  in  the  correct  locations  representing  the  joints  to  which  the  bar  is  connected. 

Note:  Since  the  bar  forces  are  computed  from  the  actual 
compatible  deflections  of  the  structure,  that  equi¬ 
librium  is  satisfied  at  each  joint  becomes  a  nec¬ 
essary  and  sufficient  condition  that  the  final  solution 
is  unique. 

Initial  Condition  of  Core:  Common  variables  IDP  and  NK  must  be  present. 


Error  Stops:  None. 

Output:  Sum  of  all  the  bar  forces  at  each  joint  are  written  on  tape  B2  for  listing. 

Note:  CHECK  is  simply  an  internal  check  of  consist¬ 
ency  of  operation.  It  indicates  that  a  structure 
has  been  analyzed  correctly;  however,  it  may 
not  be  the  desired  structure  if  some  error  has 
been  made  consistently  throughout  the  analysis, 
e.g.,  a  joint  coordinate  that  is  wrong  in  both 
STRDC  and  BARSOL  will  not  be  shown  as  an  error. 
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G.  STAIR  SERVICE  ROUTINES 


1.  Master  Routine 


(a)  Initializes  various  storage  locations  at  the  start  of  each  problem. 

(b)  Reads,  checks  and  stores  a  block  of  initial  data  for  each  problem. 

(c)  Interprets  an  input  call  for  each  step  of  an  ANALYSIS  PROGRAM; 
requests  M  ISLAM  to  read  the  routine  called  into  core,  and  transfers 
to  the  entry  point  of  that  routine. 

Error  Stops: 

(a)  NJREAC  >  125. 

(b)  10  <  NPROG  <  0  (NPROG  is  the  code  ntunber  of  a  STAIR  operating 
routine  on  a  call  card). 

(c)  Call  cards  out  of  consecutive  order  (may  be  caused  by  input  error  at 
preceding  step). 

(d)  NK  >  28. 

The  Master  Routine  serves  as  the  main  program  when  the  STAIR  System  is  loaded  by  a  BSS 
loader.  Control  passes  directly  to  it  when  transfer  vector  processing  is  completed. 

The  Master  Routine  transfer  vector  contains  the  entry  points  of  the  nine  operating  routines, 
the  MISLAM  entry  BLKIN  which  reads  processed  blocks  from  the  program  tape  and  library  I-O 
routines.  Transfer  to  each  of  the  former  routines  is  handled  by  the  standard  FORTRAN  II  CALL 
statement,  which  results  in  a  TSX  to  a  TTR  in  the  transfer  vector. 

Output:  In-core  matrices  after  any  step,  if  Sense  Switch  6  is  down. 
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2.  FIXTPE 


Operation:  Poaltlona  any  tape  on  DSC  A  to  read  a  record  m,  identified  by  a  heading  word 
32767-N  (77777.N  octal). 

Called  by:  STRK,  MATRO,  MADD,  MATIN  and  MATOUT. 

Arguments:  N,  DSC  A  logical  tape  NU,  number  of  binary  data  block  on  tape  N. 

Method:  The  first  heading  word  encountered  is  compared  to  the  desired  heading  word;  the 
tape  is  back-spaced  or  forward -spaced  according  to  the  work  input.  If  an  end-of-file  is  en¬ 
countered  in  forward-spacing  or  if  the  heading  words  are  out  of  consecutive  order,  the  tape  is 
rewound  and  forward -spaced. 

Error  Return;  II  the  tape  is  rewound  and  a  second  attempt  to  find  the  correct  heading  word 
is  not  successful,  conunon  variable  NTAPE  is  set  equal  to  one,  and  control  returns  to  the  calling 
program.  The  calling  program  senses  the  value  of  NTAPE  and  prints  an  on-line  diagnostic. 

Normal  Return:  When  a  correct  heading  word  is  read,  the  tape  is  repositioned  at  the  start 
of  the  record  containing  that  heading  word,  NTAPE  is  zeroed  and  control  returns  to  the  calling 
routine. 

Notes:  Binary  records  written  by  a  FORTRAN  program 
always  begin  with  a  zero  word,  followed  by  the 
list  in  the  "WRITE"  statement.  No  fixed-point 
variable  written  by  a  STAIR  program  may  ex¬ 
ceed  ISOOO^q;  heading  words  are  thus  safely 
distinguished  from  any  other  legal  number. 
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3.  M  A  VERT 


r 


C3peratlon: 

(a)  Expands  a  matrix  stored  in  triangular  array  at  the  head  of  STSEC  into 
square  form  with  heading  words  before  each  row  at  the  head  of  STRDC. 

(b)  Inverts  the  matrix  at  the  head  of  STRDC  (using  SHARE  routine  CLMIV-2). 

(c)  Contracts  the  inverted  matrix  back  into  triangular  array  without  heading 
words  and  stores  it  at  the  head  of  STSEC. 

Called  by:  EFFRDC  and  DEFSOL. 

Arguments:  STRDC,  STSEC,  NRA  (NRA  is  the  dimension  of  the  matrix  to  be  inverted). 

Error  Stop:  Underflow  or  overflow  during  a  floating  point  operation  causes  transfer  to  a 
floating  point  trap  routine.  If  the  number  can  not  be  normalized,  a  stop  will  occur  in  the  trap 
routine  without  an  on-line  diagnostic.  If  one  or  more  joints  among  those  to  be  EFFRDCed  or 
any  of  the  joints  in  DEFSOL  are  still  planar,  then  a  stop  will  occur  in  MAVERT  without  a 
diagnostic. 


Notes:  MAVERT  was  originally  assembled  for  IBM-704, 
and,  since  there  is  no  ii^ut  or  output,  the  binary 
deck  of  the  704  version  is  usable  on  IBM-709  or 
-7090.  Several  TQO  instructions  in  the  program 
are  treated  as  NOP's  by  the  709  and  7090,  which 
do  not  use  the  TQO  instruction. 

4.  MISLAM  for  IBM-7090 

Operation:  Processes  and  writes  a  sequence  of  routines  on  a  program  tape  from  their  re¬ 
locatable  binary  decks.  All  these  routines  must  be  in  the  form  of  standard  FORTRAN  II  sub¬ 
routines  except,  of  course,  the  main  program  (master  routine).  During  the  running  of  a  prob¬ 
lem,  certain  of  these  routines  may  be  read  back  and  stored  permanently  in  core.  Others  may 
be  read  back  individually  into  a  single  core  area  when  called  by  an  in-core  program. 

Definitions: 

(a)  Routines  stored  permanently  in  core  make  up  a  "core  block." 

(b)  Routines  stored  temporarily  in  core  are  called  "tape  blocks."  ' 

(c)  T/  cards  are  special  transfer  cards  which  separate  various  blocks 
in  the  binary  program  deck.  T/O-T/4  cards  have  a  9-row  punch 
in  column  1  and  0-4  in  the  address  field  of  the  9-left  word. 

(d)  A  "snapshot"  is  a  tape  record  which  contains  all  the  programs  written 

in  core  at  the  time  the  snapshot  is  taken.  The  lower  part  of  the  snapshot 
contains  the  MISLAM  program,  a  group  of  core-block  routines,  the  tape 
block  last  loaded  and,  possibly,  a  second  group  of  core-block  routines. 

The  upper  part  of  the  snapshot  contains  the  BSS  loader  and  program 
symbol  table. 

(e)  A  "bootstrap"  is  a  self-loading  program  written  as  the  first  record  on 
the  program  tape.  The  lower  part  of  the  bootstrap  locates  the  last  snap¬ 
shot  written  on  the  program  tape.  The  upper  part  reads  that  snapshot 
into  core  and  transfers  to  the  entry  point  ot  the  MISLAM  program. 

(f )  "Dummy  records"  are  two-word  records  separating  bootstrap,  snapshot 
ancf  tape -block  records.  Their  presence  allows  some  saving  in  time 
during  backspacing.  Two  words  are  needed  to  allow  a  program  delay 
during  a  tape  search. 

(g)  An  "end  record"  is  a  two-word  record  which  marks  the  last  information 
written  on  the  program  tape.  No  reading  operation  may  pass  an  end  record. 
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Functioning  of  T/  Cards:  When  the  MISLAM  program  or  a  anapahot  la  read  Into  core,  con> 
trol  paaaea  almoat  at  once  back  to  the  BSS*  loader.  The  BSS  reada  cards  from  the  on-line 
card  reader.  Instruction  and  program  cards  are  Interpreted  In  the  normal  manner  with  the  cur¬ 
rent  relocation.  When  a  transfer  card  is  read,  control  passes  back  to  MISLAM,  which  branches 
to  one  of  the  following  options: 

a.  T/1  Card 

A  T/l  card  initiates  a  tape  block. 

If  the  preceding  block  was  a  tape  block  {i.e.,  if  another  T/l  card  was  read  previously),  that 
block  is  written  on  tape  with  its  transfer  vector  unprocessed.  Relocation  is  reset  to  the  begin¬ 
ning  of  the  tape-block  area,  the  tape-block  area  is  cleared  and  control  returns  to  the  BSS. 

If  the  preceding  block  was  not  a  tape  block  (it  must  then  be  the  lower  core  block),  the  start¬ 
ing  location  of  all  tape  blocks  is  set  before  control  passes  back  to  the  BSS. 

The  length  of  the  tape-block  area  is  set  with  the  length  of  the  longest  tape  block  read. 

b.  T/2  Card 

A  T/2  card  initiates  a  core  block. 

If  the  preceding  block  was  a  tape  block,  it  is  written  unprocessed  on  tape.  Relocation  is 
set  to  the  end  of  the  longest  tape  block,  and  control  returns  to  the  BSS. 

c.  T/3  card 

A  snapshot  is  written  on  the  program  tape,  and  control  returns  to  the  BSS. 

d.  T/0  card 

This  is  the  normal  FORTRAN  II  transfer  card. 

If  the  preceding  block  was  a  tape  block  (i.e.,  if  the  preceding  T/  card  was  a  T/l  or  a  T/2 
preceded  by  a  T/l),  it  is  written  unprocessed  on  the  program  tape. 

The  program  tape  is  rewound  and  positioned  after  the  bootstrap  record.  Successive  tape 
blocks  are  read  from  the  tape,  their  transfer  vectors  are  processed  and  they  are  rewritten  as 
processed  tape  blocks  beyond  the  present  end  record.  Core  block  programs  are  similarly 
processed.  Finally,  the  symbol  table  and  BSS  are  cleared,  and  control  transfers  to  the  main 
program  (master  routine). 

e.  T/4  Card 

If  all  tape  blocks  are  correctly  processed  and  written  on  tape,  only  the  core  block  is  proc¬ 
essed  before  control  passes  to  the  main  program.  Usually  this  would  be  used  on  all  normal 
runs  instead  of  a  T/O  card,  except  those  which  involve  writing  of  tape  blocks. 


♦  Note:  The  BSS  loader  for  the  709  or  7090  must  be  a  10-card  loader  due  to  timing  in  the  in¬ 
structions. 
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ntOGRAM  TAPE  FORMATi 


BEGINNING  OF  TAPE  BOOTSTRAP 

RECORD  GAP  (TYPICAL) 

DUMMY  RECORD 

UNPROCESSED  TAPE  BLOCK  I 

DUMMY  RECORD 

UNPROCESSED  TAPE  BLOCK  2 

DUMMY  RECORD 

UNPROCESSED  TAPE  BLOCKS 

DUMMY  RECORD 

SNAPSHOT 

DUMMY  RECORD 

UNPROCESSED  TAPE  BLOCK  N 
DUMMY  RECORD 

PROCESSED  TAPE  BLOCK  1 

DUMMY  RECORD 

PROCESSED  TAPE  BLOCKS 

DUMMY  RECORD 

PROCESSED  TAPE  BLOCK  N 

DUMMY  RECORD 

END  RECORD 


RECORD  FORMATS: 


A.  UnprocMMd  top*  block 


CODE  word  (1000  +  block  numbor) 
dummy  word 

unprocoMod  routino 

CODE  word  (2000  +  block  numbor) 
procoued  routino 


C.  Snopihot 


mopihot  codo  word  (1000  mopthot  No.) 
quantity  of  wordt  in  lowor  port  of  oiopiiiot 

lowor  part  of  anapthot 

rood  command  for  uppor  port  of  mqpihot 
dummy  word 


uppor  port  of  arapihot 


D.  Dummy  lUeord 


77776  octal 
dummy  word 


77777  octal 
dummy  word 


Notai:  tdonrtfying  code  numbon  oro  hold  In  Iho  docromont  of  tho  fint  word  of 
each  record.  Dummy  wcrdi  oro  required  for  proQrom  deloy  while  the 
preceding  vrord  It  being  rood. 
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Reading  Back  Tape  Blocks:  Processed  tape  blocks  are  read  Into  core  by  MISLAM  at  its 
BLKIN  entry  point.  This  entry  ntay  be  made  from  any  core-block  program  with  a  CALL 
BLKIN  (I)  statement.  The  tape  is  positioned  to  read  tape  block  I,  the  block  is  read  and  control 
passes  back  to  the  calling  program.  The  tape  block  is  entered  by  a  second  CALL  statement. 

Assembly  of  MISLAM  for  IBM-709  or  7090:  The  SAP-coded  version  of  MISLAM  for  704  re¬ 
quires  several  modifications  in  the  assembly  program  which  allow  various  parts  of  MISLAM  to 
overlap.  (In  particular,  after  the  bootstrap  record  is  written,  the  main  part  of  MISLAM  is 
stored  over  the  bootstrap  in  core.)  These  modifications  could  easily  be  made  in  9AP,  so  certain 
709  and  7090  operation  codes  are  defined  in  SAP  language  and  the  program  is  assembled  on  704, 
using  CORSAP  as  Assembler. 

Error  Stops: 

(a)  Tape  failures  during  reading  or  writing  operations. 

(b)  No  snapshot  written  when  writing  program  tape. 

(c)  Proper  sequence  of  T/  cards  violated.  Proper  sequence  when  writing 
program  tape  is  T/l  —  T/2  —  T/O.  T/3  cards  may  be  given  at  any  time; 
a  t74  card  may  only  be  used  alone. 

(d)  Number  of  tape  blocks  >  20. 

(e)  Nvimber  of  snapshots  >1000. 

(f )  Reference  from  one  tape  block  to  another  (any  program  in  the  core  block 
may  refer  to  a  tape  block  when  it  is  read  in). 

Writing  the  STAIR  Program  Tape:  The  FORTRAN  II  and  SCAT  routines  comprising  the 
STAIR  system  are  compiled  into  relocatable  binary  decks  in  the  usual  manner.  These  decks 
are  combined  with  T/  cards  into  the  following  deck  which  Is  read  on-line  and  written  on  tape  B4. 
The  order  of  operating  routines  must  be  followed  exactly  as  shown.  Since  control  passes  to  the 
main  program  following  the  T/O  card,  data  may  be  waiting  for  an  analysis  program  when  the 
program  tape  is  written. 


klnoiy  Oaek  S«qu«ne« 
fer  STAR  Program  Top* 


[=3 


FORTRAN  BSS  LOADER* 

MISLAM 

MAIN 

LIBRARY  SERVICE  ROUTINES 

STAIR  ROUTINES  FIXTPE,  MAVERT 
T/1  CARD 

STAIR  ROUTINE  STRIX 
T/1  CARD 

STAIR  ROUTINE  EFFRIX 
T/1  CARD 

STAIR  ROUTINE  MATRO 
T/1  CARD 

STAIR  ROUTINE  MADD 


T/1  CARD 

STAIR  ROUTINE  MATOUT 
T/1  CARD 

STAIR  ROUTINE  MATIN 
T/1  CARD 

STAIR  ROUTINE  DEFSOL 

T/1  CARD 

STAIR  ROUTINE  BARSOL 


T/1  CARD 

STAIR  ROUTINE  CHECK 

T/3  CARD 
T/0  CARD 
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Subaequent  Runs  of  Program:  When  a  program  is  written  on  tape  by  MISLAM,  there  are 
several  ways  of  initiating  a  run  of  the  program. 

(a)  New  program  is  to  be  added  to  core  block.  (A  new  version  of  a  program 
will  rep^ce  a  prior  version  of  the  program  in  the  transfer  vectors  of  all 
calling  programs.) 

-  get  tape  B4  card  (loads  bootstrap) 

-  T/2  card 

new  core-block  program 

-  T/3  card  (optional) 

-  T/O  card  ^ 

(b)  No  new  program  to  be  read;  tape-block  transfer  vectors  are  already  processed. 
(This  is  the  normal  procedure  for  initiating  a  run  when  the  program  is  debugged.) 

-  get  tape  B4  card 

-  T/4  card 

H.  STAIR  LOAD  PROGRAMS 

Operation:  Form  load  matrices  for  each  unit  of  a  structure  and  write  these  matrices  as 
consecutive  binary  records  on  tape,  which  may  be  read  back  by  the  STAIR  program  during  an 
analysis  problem. 


1.  SLOP  1 

Method:  For  each  of  a  number  of  loading  conditions,  joint  loads  are  read  into  core  and 
stored  (in  a  format  similar  to  displacement  storage  by  DEFSOL).  When  all  loads  for  all  load¬ 
ings  are  stored,  joint  lists  identifying  the  matrices  of  consecutive  units  are  read,  and  load 
matrices  are  assembled  according  to  the  order  of  joints  in  the  joint  list.  Rows  of  each  load 
matrix  correspond  to  forbe  components  and  columns  to  loading  conditions.  The  matrices  are 
stored  columnwise.  An  end-of-file  is  written  after  the  last  load  matrix.  A  heading  word  is 
written  as  the  first  word  in  each  load  matrix  record  (32767-u,  where  u  is  the  number  of  a  unit). 
As  each  term  is  taken  from  the  load  storage  area  to  be  stored  in  the  matrix  assembly  area,  its 
register  in  the  load  storage  area  is  cleared;  this  prevents  duplication  of  loads  at  joints  shared 
by  two  or  more  units. 

Error  Stops: 

(a)  3  X  number  of  loadings  x  number  of  joints  in  structure  >  26000. 

(b)  3  X  number  of  loadings  x  number  of  joints  in  any  unit  >  2000. 

(c)  Joint  lists  out  of  consecutive  order. 

Error  Diagnostic,  but  no  Stop: 

(a)  Illegal  joint  number  with  input  loads. 

(b)  Illegal  joint  number  in  a  unit  joint  list. 

Output;  Load  matrices  are  always  stored  as  binary  records  on  tape  AS.  If  Sense  Switch  3 
is  depressed,  load  matrices  will  also  be  written  on  tape  B2  in  form  for  listing. 
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storage  Addreaalng:  Joint  numbers  have  the  form  1  j,  where  i  Is  a  unit  number  and  j  is  a 
joint  number  within  that  unit.  Addresses  in  the  load  storage  area  are  computed  with  a  "unit 

i  j 

addressing  array"  lUADR,  where  lUADR^  =  £  3  Uj^  (Uj^  is  the  munber  of  joint  numbers  pre¬ 
fixed  with  unit  number  k) .  The  location  A  of  the  x-component  of  load  at  joint  i  j  for  loading  1 
is  thus:  A  =  (1  -  1)  N  +  lUADR^  4  3(j  -  1)  4  1(N  is  3  x  number  of  joints  in  the  structure).  An 
illegal  joint  number  is  either  one  of  the  following; 

(a)  i  >  total  number  of  units. 

(b)  lUADRj  4  3(j  -  1)  4  1  >  lUADRj  4  1. 
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2.  SLOP  2 

Method;  In  addition  to  the  functions  performed  by  SLOP  1,  SLOP  2  computes  joint  loads 
due  to  dead  weight  of  the  structure.  Input  to  the  program  is  divided  into  groups;  the  first  word 
in  each  group  selects  the  function  to  be  performed  (store  hand-computed  loads,  compute  dead 
loads  or  assemble  load  matrices). 

Optional  Error  Stops:  (no  stop  if  SSI  is  down) 

(a)  Illegal  joint  numbers  in  input  list  of  coordinates. 

(b)  Illegal  bar  number  in  input  list. 

(c)  Illegal  joint  number  in  input  list  of  loads. 

(d)  Illegal  joint  number  in  the  joint  list  of  a  unit. 

Unconditional  Error  Stops: 

(a)  Illegal  input  code  at  head  of  an  input  group. 

(b)  Illegal  loading  number  on  input  group  heading  card. 

(c)  3  X  niunber  of  loadings  x  number  of  joints  in  structure  >  19000. 

(d)  Unit  joint  lists  out  of  consecutive  order. 

(e)  3  X  number  of  loadings  x  number  of  joints  in  any  unit  >  2000. 

Output:  In  addition  to  the  output  options  of  SLOP  1,  total  weight  of  the  structure  and  co¬ 
ordinates  of  the  center  of  gravity  are  printed  out  on-line  for  each  dead  load  calculation. 
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CHAPTER  5 

LIMITATIONS  AND  EXTENSIONS 


At  this  writing,  the  lull  capacities  of  the  STAIR  program  have  not  been  realized.  Theoreti¬ 
cally,  a  structure  that  has  a  total  of  3875  joints  (including  125  reaction  joints)  can  be  analyzed 
for  one  loading  condition,  but  the  largest  structure  analyzed  tb  date  has  only  600  joints  and  three 
loading  conditions. 

For  problems  involving  large  structures  (more  than  400-500  joints),  the  limiting  restrictions 
of  the  various  routines  can  be  circumvented  as  follows: 


SLOP: 


STRDC: 

EFFRIX: 

MATRCh 

MADD: 


MATOUT: 

MATIN: 

DEFSOL: 

BARSOL: 


If  there  are  no  gravity  loads,  SLOP  represents  no  restriction. 

If  gravity  loads  are  wanted,  the  structure  can  be  bro)cen  into 
several  sections  and  the  gravity  loads  computed  for  each  sec¬ 
tion  separately.  Then,  by  superimposing  the  loads,  SLOP  can 
be  run  again  with  the  gravity  loads  as  applied  loads.  Note  that 
storage  locations  also  could  be  reassigned  in  SLOP  to  accom¬ 
modate  a  larger  problem,  but  this  measure  is  somewhat  limited. 

The  structure  can  be  brolcen  into  any  number  of  units. 

This  routine  can  be  called  any  number  of  times  if  necessary.* 
This  is  limited  by  STRDC,  and  thus  is  not  a  limitation  in  itself. 

EFFRIX  can  be  called  any  number  of  times  before  a  MADD  op¬ 
eration;  thus,  the  matrices  to  l>e  added  can  be  reduced  to  the 
prqier  size. 

This  is  limited  only  by  the  limitations  of  the  other  subprograms. 
This  is  limited  only  by  MATOUT. 

EFFRIX  can  be  called  any  number  of  times  before  a  DEFSOL 
operation.  The  only  other  limitation  is  discussed  below. 

This  routine  can  be  called  any  number  of  times. 


CHECK:  The  restrictions  cf  this  routine  do  not  affect  the  program. 


The  fourth  restriction  (see  p.  39)  of  DEFSOL  represents  the  only  limitation  on  the  STAIR 
program.  If  this  restriction  becomes  too  severe,  the  program  could  be  rewritten  to  accommo¬ 
date  a  special  problem,  but  an  experienced  programmer  should  be  consulted. 

It  should  be  Icept  in  mind  that  the  STAIR  program  solves  only  linear  problems  since  the 
stiffness  matrix  of  the  structure  is  based  on  the  undeformed  geometry.  If  deflections  are  so 
large  that  the  angles  of  connection  of  members  at  joints  are  altered,  the  results  should  be  treated 
as  only  approximations  of  the  actual  structural  behavior  since  the  problem  is  no  longer  linear. 
Further,  there  is  no  provision  in  the  program  to  handle  budding  of  individual  members  of  the 
structure  or  buckling  of  the  structure  as  a  whole.  Consequently,  structural  instability  other  than 
geometric  instability  must  be  analyzed  by  other  methods,  i.e.,  buckling  of  a  trussed  arch  or 
column  would  not  be  indicated  in  the  results  of  a  STAIR  ANALYSIS  PRCXIRAM. 

Computer  round-off  error  becomes  important  when  solving  extremely  large  problems.  It 
is  difficult  to  say  how  much  round-off  error  will  affect  any  one  problem,  since  the  parameters 
affecting  the  error  are  a  function  of  the  problem  to  be  solved.  In  general,  the  types  of  structures 
analyzed  hy  STAIR  are  not  greatly  affected  by  round-off  error,  but  it  cannot  be  discounted  as  a 
possible  source  of  inaccuracy  in  the  results.  For  any  given  problem,  it  is  possible  to  evaluate 
the  round-off  error  if  an  experienced  programmer  or  mathematician  is  consulted.  Experience 


*EFFR1X  can  be  called  only  cmce  for  a  given  unit  if  MATRO  later  forms  another  unit  from  the  first. 
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in  problem*  with  up  to  200  joint*  has  shown  the  results  were  good  to  five  iiiTi.’i'/can*  mres 
However,  this  should  not  be  used  as  a  standard  of  ctmparison,  but  only  as  r  | 

In  its  present  state,  the  program  will  solve  for  Reflections  for  only  gri.vity  and  i 
and  the  structuire  has  to  be  made  of  rtietnbers  of  constant  density  and  moduius  of  elasticity^ 
ever,  deflections  caused  by  temperatur';  changes  or  fa^ication  errors  can  >e  handled  with  i 
modifications  of  the  program.  By  making  appropriate  chaises  in  the  areas  of  bars  in  SLOP,  or 
the  actual  STAIR  ANALYSIS  PROGRAV.:.  members  whichfhave  different  densities  or  moduli  of 
elasticity  can  be  handled  by  the  compute  i-.  ^  ^  t 

Occasionally,  reactions  may  be  in  li  rectionr  i.^ther  than  along  the  cooc^nate  axes, 
problem  can  be  easily  handled  by  us:.ng  i  simple  o^ie-  or  twp-bar  "tripod"  with  short  members 
of  large  area.  The  bar(s)  of  the  tripod  ire  oriented  in  the  (direction  of  thi'  desired  reaction  and  ^ 
the  standard  x,  y  and  z  reactions  iiold  the  ends  of  the  ba^s  opposite  the  structure.  '  i  % 
Other  changes  in  the  input  data  or  t  ie  program  itself  ^ill  extend  the  .'NALYSIS  PROGRAM 
to  cover  any  type  of  loading  or  structurc.l  action  to  which  a  pin-jointed  st  r  icture  may  be  subjected. 
In  making  modifications  to  the  input  datr  >r  the  program  itself,  jk  clear  i  i  derstanding  of  the 
mathematical  method  of  solution  bl  the  lil’AIR  program  is  necessary.  Although  some  analogies 
exist  between  the  matrix  operations  invoh'ed  in  the  solution  of  a  problem  and  physical  operations 
In  the  actual  structure,  these  analogies  sometifffes  break  down  and  will  c  :use  erroneous  results. 
The  mathematical  model  of  the  structure  ]nay  not  always  b’h^n  exact  rep:  oduction  of  the  physical 
structure,  and  it  cannot  be  overemphasized  that  sulitle  differ^ces  in  the  liehavior  of  the  struejurg 
and  its  mathematical  model  may  cause  baffling  and  erroneous  fsosults.  •" 
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